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ABSTRACT 


The United States Navy (USN) uses Optimal Track Ship Routing provided by ship 
routing officers (SRO) to aid in the safe transit of its ships. When a ship makes a transit, 
the ship provides the SRO an origin, a destination, and a date of departure, and the SRO 
will generate a route for the ship to proceed along. Avoiding severe weather is the most 
important consideration in determining the route. In addition to safe transit, the USN also 
focuses on fuel efficiency. In recent years, the meteorology and oceanography 
community has been providing more products that estimate the uncertainty in 
environmental forecasts. However, it is not known how much that uncertainty affects or 
should affect ship routing. This thesis explores the sensitivity and robustness of 
optimized ship routes generated by the Ship Track and Routing System optimizer to 
uncertainty in the environment. 
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EXECUTIVE SUMMARY 


Meteorology and oceanography (METOC) forecasters rely heavily on numerical models. 
Numerical weather prediction (NWP) models are discretized mathematical models 
consisting of basic differential equations simulating the fluid dynamic process and 
physics of the environment. Based on Newton’s second law of motion, basic equations 
are solved for any point (location and altitude) at any given time in the atmosphere for 
dependent variables that describe the atmosphere and stepped through time to obtain 
predictions of the future behavior of the atmosphere. Depending on the horizontal and 
vertical spacing of the discretization, some aspects of the equations are solved explicitly 
while smaller-scale operations are “parameterized.” Solutions to all equations for each 
point are coupled with neighboring points, and all equations with their solutions are 
stepped forward in time simultaneously until the desired forecast length is achieved. 

Ensembles are a way of using numerical models to quantify uncertainty in 
environmental conditions. An ensemble forecast can be generated from the use of two or 
more different forecast techniques or numerical models, or using multiple numerical 
models with different initial conditions. 

NWP forecasts are driven primarily by initial conditions. Initial conditions are 
the model’s variable values at t=t 0 , where t represents time, from which the model’s 
forecast is stepped forward. The initial estimates of the model variables are more 
important than the model boundary conditions. The initial value of every variable at 
every location in the model must be specified. Since observations are not available in all 
locations and because the available observations often differ from the variables solved by 
the NWP model, initial conditions are uncertain. To account for this uncertainty, 
ensemble forecasting runs the model several times. Each time the model is run, the initial 
conditions are perturbed slightly. Small changes in initial conditions can lead to large 
differences in forecast values, and the collection of ensemble forecasts provides an 
estimate of the expected error in the forecast. 
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This thesis explores the sensitivity and robustness of optimized ship routes 
generated by the Ship Track and Routing System (STARS) optimizer to uncertainty in 
the enviromnent. Environmental uncertainty is represented by an ensemble generated 
from numerical weather models. STARS uses the Naval Operational Global Atmospheric 
Prediction System (NOGAPS) for the wind model. The Wavewatch III (WWIII) model 
is used as the wave model, and the global Navy Coastal Ocean Model (NCOM) is used as 
the current model. 

This thesis sets out to answer questions regarding the sensitivity and robustness of 
ship routes to environmental conditions and how important it is to take into consideration 
environmental uncertainty when determining a robust ship route. The sensitivity of the 
optimized routes showed that there were few routes that differed greatly in distance 
varying environmental forecasted conditions. The robustness of the optimized routes 
showed that the routes were safe under various environmental conditions. The time of 
year, with respect to seasons, is a major factor in the uncertainty of environmental 
conditions. The time of year and the forecasted weather conditions should definitely be 
taken into consideration when trying to determine a robust ship route. 

Given that weather is a variable that cannot be controlled, it is important to 
understand how that uncertainty can affect a ship route. The optimized routes that were 
generated by STARS showed little sensitivity and were robust to the environmental 
conditions. Although the environmental conditions proved to have an impact on the 
optimized routes, the generated routes were, for the most part, the best and the safest 
route to take. However, because forecasted weather is not a guarantee but an estimate of 
a possible occurrence, ship routes based on forecasted environmental conditions are 
subject to change with the changing conditions. The optimized routes can be used as a 
basis, but it should be clear that even an optimized route may need to be altered to 
maintain robustness. 
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I. INTRODUCTION 


A. BACKGROUND 

The mission of the United States Navy (USN) is to provide a naval force capable 
of ensuring that the freedom of the seas is maintained for the United States and its allies. 
Moored pier side or underway, ships are vulnerable to the environmental conditions that 
they encounter. A simple transit across the Atlantic Ocean can easily become a rough 
voyage if the ship encounters high winds, which in turn will cause a high sea state. 

Environmental conditions play a crucial role in safe navigation. Each USN ship 
class has a set of safe operating envelopes (SOE). These SOE indicate the maximum seas 
and the maximum winds in which the ship can safely operate. If conditions that exceed 
the SOE are not avoided, the ship is at risk to receive extensive damage or even be a total 
loss; there is also a risk of injury to or loss of ship crew members. 

History has shown that extreme environmental conditions can cause damage to a 
ship and the crew that it holds. The Naval Historical Center lists approximately forty 
weather-related incidents in naval history from 1781 to 1949, most of which were caused 
by gale force winds, heavy squalls, hurricanes, stonns, and typhoons (Naval Historical 
Center, 2005). In December 1944, the ships of Task Force 38 were overwhelmed by a 
typhoon that claimed three destroyers and the lives of approximately 790 officers and 
men; the typhoon also caused severe damage to a cruiser, five aircraft carriers and their 
planes, and three other destroyers, and injured over eighty personnel (Naval Historical 
Center, 2001). More recently, in November 2005 while returning home from a UNITAS 
deployment, the USS Ross (DDG51) was hit by rough seas, causing a crack in the hull 
and damaging the ship’s sonar equipment rooms (A. Corbin, personal communication, 
April 13, 2011). In each of these incidents, the effects of the unexpected environmental 
conditions encountered proved costly. 

The USN uses Optimal Track Ship Routing (OTSR) provided by ship routing 
officers to aid in the safe transit of its ships. Typically, ship routing officers (SRO) 
monitor as many as 5,300 routes a year (F. Sullivan, personal communication, April 15, 
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2011). Prior to departure, a ship provides the SRO an origin, a destination, and a date of 
departure, and the SRO will generate a route for the ship to proceed along that is 
forecasted to be consistent with the ship’s operational limits. However, the route is 
subject to change. 

There are many reasons for a ship to change its route, but the most typical 
weather-related reason to alter the route is to avoid or minimize its exposure to sea and 
wind conditions which exceed the ship’s SOE. Even though the ship receives weather 
advisories, there are some cases in which the advisories are not sent to or received by the 
ship in time to prepare for the change in the weather. Uncertainties in weather play a 
vital role in the safe navigation of the ship. The SRO provides an alternative route to 
modify a ship’s route due to hazardous weather or tropical cyclone avoidance. If a USN 
ship chooses to be non-compliant with the OTSR diversion recommendations, they must 
receive permission from the respective Fleet Commander via the ship’s respective fleet 
number. This situation usually occurs when the mission is of great importance. 

A second important consideration in ship routing is fuel. The Ticonderoga-class 
(CG47) cruiser has the capacity to hold over 15,000 barrels (bbls) of distillate marine fuel 
(DFM) (United States Department of the Navy [USDoN], 2007). As stated in USDoN 
(2007), the petroleum, oil and lubricant (POL) requirement for the CG47 class is over 
750 bbls per day for sustainment at sea. This is an average POL fuel consumption rate, 
but there are other factors that influence fuel consumption. Fuel consumption is affected 
by speed traveled, distance traveled, and the environmental conditions experienced along 
the route. When considering factors that influence fuel consumption, the intensity and 
direction of both waves and winds are primary (Brown, Kline, Rosenthal & Washburn, 
2007). 

Meteorology and oceanography (METOC) forecasters rely heavily on numerical 

models. Numerical weather prediction (NWP) models are discretized mathematical 

models consisting of basic differential equations simulating the fluid dynamic process 

and physics of the environment. Based on Newton’s second law of motion, basic 

equations are solved for any point (location and altitude) at any given time in the 

atmosphere for dependent variables that describe the atmosphere and stepped through 

2 



time to obtain predictions of the future behavior of the atmosphere (Lynch, 2008). 
Depending on the horizontal and vertical spacing of the discretization, some aspects of 
the equations are solved explicitly while smaller-scale operations are “parameterized.” 
Solutions to all equations for each point are coupled with neighboring points, and all 
equations with their solutions are stepped forward in time simultaneously until the 
desired forecast length is achieved. 

Ensembles are a way of using numerical models to quantify uncertainty in 
environmental conditions. Sivillo, Ahlquist, and Toth (1997) define an ensemble forecast 
as “a collection of two or more forecasts that verify at the same time.” An ensemble 
forecast can be generated from the use of two or more different forecast techniques or 
numerical models, or using multiple numerical models with different initial conditions. 

NWP forecasts are driven primarily by initial conditions. Initial conditions are 
the model’s variable values at / =/ 0 , where t represents time, from which the model’s 
forecast is stepped forward. The initial estimate of the model variables is more important 
than the model boundary conditions. The initial value of every variable at every location 
in the model must be specified. Since observations are not available in all locations and 
because the available observations often differ from the variables solved by the NWP 
model, initial conditions are uncertain. To account for this uncertainty, ensemble 
forecasting runs the model several times. Each time the model is run, the initial 
conditions are perturbed slightly. Small changes in initial conditions can lead to large 
differences in forecast values, and the collection of ensemble forecasts provides an 
estimate of the expected error in the forecast. 

B. RESEARCH QUESTIONS 

An important question for the METOC community is how much METOC 
forecasts affect ship routes and how severe are the consequences of following a non- 
optimal route. Therefore, this thesis addresses the following questions: 

1. How sensitive are optimal ship routes to environmental conditions, i.e., 
waves and wind as predicted by distinct forecasts? 
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2. How robust is the optimized ship route to the variability of environmental 
conditions as represented by multiple distinct forecasts for the same time 
period? 

3. How valuable is taking environmental uncertainty into account when 
finding a robust ship route? 

C. LITERATURE REVIEW 

There has been little research done on the effect of uncertainty in environmental 
conditions on naval operations. Stoughton (2010) focuses on the wind thresholds as they 
apply to a USN nuclear powered aircraft carrier (CVN) during ammunition offload. 
Stoughton uses an ensemble to generate a probabilistic forecast and bases a routing 
decision on that forecast. In this thesis, each ensemble member is treated as a separate 
deterministic forecast and an optimal route is found with respect to each of those 
forecasts. 

Other authors have addressed optimizing ship routes with respect to the 
environment. Montes (2005) uses shortest path models to analyze ship routing and 
provide an operations research tool to aid in the initial route planning, recommendation, 
and diversions. Montes shows that the use of OTSR is a useful tool to analyze parameter 
sensitivities when planning route diversions. Vlachos (2004) proposes two methods for 
generating nearly optimal ship routes, using transit time and safety (a function of seas and 
winds) as the objective. Vlachos treats the environment as deterministic. Dolinskaya, 
Kotinis, Parsons, and Smith (2009) investigate the involuntary loss of speed caused by 
added resistance and the limitations imposed on the vessel by operational constraints. 
The results of the seakeeping model provide the maximum attainable speed for a given 
sea state and ship heading. 
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D. 


SCOPE OF THESIS 


This thesis explores the sensitivity and robustness of optimized ship routes 
generated by the Ship Track and Routing System (STARS) optimizer to uncertainty in 
the enviromnent. Environmental uncertainty is represented by an ensemble generated 
from numerical weather models. STARS uses the Naval Operational Global 
Atmospheric Prediction System (NOGAPS) for the wind model. The Wavewatch III 
(WWIII) model is used as the wave model, and the global Navy Coastal Ocean Model 
(NCOM) is used as the current model. 
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II. DESIGN OF EXPERIMENTS 


Using optimization code provided by the Naval Research Laboratory (NRL) with 
STARS, “optimal” ship routes were generated for multiple forecast environments for 
three frequently traveled origin-destination pairs from June 2010 to December 2010. 
These ship routes were generated to explore the sensitivity and robustness of optimal 
routes to uncertain environmental conditions as modeled with ensemble wind, wave, and 
current forecasts. 

A. SIMULATION ENVIRONMENT 

STARS was developed by Oceans System Incorporated in Alameda, California. 
It is a software package that generates optimal ship routes based on weather model data 
(Montes, 2005). It takes as an input forecast METOC conditions and produces a route 
that minimizes total energy required for propulsion subject to constraints on METOC 
environmental conditions for a given forecast environment (J. Etro, personal 
communication, May 13, 2011). Appendix A contains a sample of a STARS input file, 
which shows the environmental constraints applied to the generated routes. These 
constraints are 35 knots of wind speed, 65 knots of relative wind speed, and 25-foot seas. 
The STARS program uses the DD963 Spruance-class destroyer as its modeled platform. 
Although this class of ship has been decommissioned, its characteristics have been 
maintained in today’s Ticonderoga-class cruiser. When STARS was presented to the 
USN, it was part of the Automated Optimum Track Ship Routing (AOTSR) system, a 
system used by SRO to assist them in OTSR. However, ship routers were reluctant to use 
this software because it was designed for merchant vessels and was difficult to use 
(Montes, 2005). 

1. METOC Environments 

The environments in STARS used in this thesis are the members of a sixteen- 

member ensemble generated from the NOGAPS model for winds, WWIII model for 

waves, and NCOM for currents. The initial conditions for the sixteen environments are 

generated by perturbing the most recent analysis using an ensemble transform analysis 

7 



perturbation scheme (McLay, Bishop, & Reynolds, 2010). Hogan, Rosmond, and Gelaro 
(1991) give a detailed description of the NOGAPS model. The NOGAPS model is 
crucial to nearly every environmental product produced by the Fleet Numerical 
Meteorology and Oceanography Center (FNMOC), and it is the foundation of the Navy’s 
operational and environmental central site fleet support. In STARS, WWIII is an offline 
model. It can receive any wind field forecast, which is then used to detennine waves. 
The wind fields are generated in the NOGAPS model and passed to the WWIII model. 
The global NCOM provides surface boundary conditions, which in turn are provided by 
the NOGAPS (National Weather Service [NWS], 2009). NCOM uses moisture and heat 
from the lower atmosphere to generate currents and receives atmospheric forces from the 
NOGAPS (J. Hansen, personal communication, May 16, 2011). The WWIII waves and 
NCOM currents are generated for each NOGAPS environment, providing a forecast 
environment that includes winds, waves, and currents for each of the sixteen 
environments. 

The environments are arbitrarily named and indexed i = = 16. All are 

generated the same way. There are two other environment options available in STARS. 
The “analysis” is the best estimate of the environmental; it is also the time zero forecast 
for the numerical model. The “ensembleaverage” environment is the average value of all 
the selected raw ensemble members (J. Hansen, personal communication, May 12, 2011). 

2. STARS Input 

STARS optimizes on a grid that is contained within a user-defined envelope. For 
each route, the input file specifies the upper bounds (UBs) and lower bounds (LBs) of the 
envelope, each represented in the form of a latitude-longitude pair to the tenth of a 
degree. UB NUMBER 01 and LB NUMBER 01 are the same latitude-longitude pair, 
which is the route origin, while UB NUMBER 04 and LB NUMBER 04 are the same and 
represent the route destination. In addition, STARS requires two UBs and two LBs that 
create the envelope. The envelope should resemble the shape of a football and it should 
not include any land masses. The envelope for each origin-destination pair was chosen 
manually. The UBs and LBs chosen were far away from the expected route as possible 
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so as to not constrain the optimized route. Figures 1 and 2 illustrate the envelope made 
by the UBs and LBs. Figure 1 is in open ocean, and Figure 2 shows how the envelope is 
adjusted to avoid land masses within the envelope. 



Figure 1. UBs and LBs Envelope for a Route Between San Francisco, CA and Pearl 

Harbor, HI (From Google Earth, 2010) 



Figure 2. UBs and LBs Envelope for a Route Between Capetown, South Africa and 

Hobart, Australia (From Google Earth, 2010) 


9 











Appendix B contains an example of the executable file for STARS. The data in 
this file are used to run each simulation, and it is this file that points to the applicable 
input file. The departure date and arrival date are data input fields that are indicated by a 
dollar sign ($) on the input file (Appendix A). The most important data for the execution 
file are the Routex environment and Weax environment data fields. Within the 
designated Routex environment, STARS takes the latitude-longitude pairs for the origin 
and destination, creates a grid-like diagram, and then optimizes the route for the given 
time of departure and arrival and the Routex environment. 

3. STARS Output 

Routes (denoted/-) are uniquely determined by origin-destination pair, start date 
(arrival date is not a constraint), and Routex environment (indexed/). A route r consists 
of a set of N r waypoints, w n ,n = \,...,N r , where the n th waypoint is defined by the triplet 

w n = (d n ,x n ,y n ) (Table 1), where: 


d n 

date-time group of n' h waypoint 

X n 

latitude of n' h waypoint, -90.0° to 90.0°, where + indicates 


the Northern Hemisphere and - indicates the Southern 


Hemisphere. 

y n 

longitude of waypoint, 0.0° to 360.0° 


Table 1. Notation and Definitions for Geographical Locations 


The notations*(z),x* (/), y„ (/) will be used to denote the STARS-generated optimal 

route latitude and longitude for the n"‘ waypoint, optimized relative to environment /. 
The STARS output for an optimized route also includes the following (Table 2): 


S n 

ship speed at n th waypoint in knots 

C n 

ship course at n th waypoint, 0.0° to 360.0° 

dist n 

distance in nautical miles (nm) to next waypoint 


Table 2. Notation and Definitions for Ship’s Course, Speed and Distance 
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Appendix C contains an example of a HyperText Markup Language (HTML) output file 
for STARS. Each row of this table corresponds to a waypoint, which are usually 
separated by six hours. 

STARS also evaluates the route with respect to an environment—called the Weax 
environment—and indexed j that does not necessarily coincide with the Routex 
environment. STARS reports the following (Table 3): 


wind _ speed ( w n , j ) 

Surface wind speed in knots for the n lh waypoint 
evaluated with respect to environment j 

wind _ direction ( w n , j) 

Surface wind direction, 0.0° to 360.0° 

sig _ wave _ height ( w n , j) 

Sea surface significant wave height in feet 

sea _ height (w n , j) 

Sea height in feet 

sea _ period ( w n , j) 

Sea period in seconds 

sea_direction (w n , j) 

Sea direction, 0.0° to 360.0° 

swel l _ heigh t(w n ,j ) 

Swell height in feet 

swell period (w n ,y) 

Swell period in seconds 

swell _ direction ( w n , j ) 

Swell direction, 0.0° to 360.0° 

current _ speed ( w n , j ) 

Current speed in feet/second 

current _ direction ( w n , j ) 

Current direction, 0.0° to 360.0° 

horsepower ( w n , j) 

HP in kilo horsepower-hours since last waypoint 


Table 3. Notations and Definitions of STARS Outputs 
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B. 


EXPERIMENTAL ROUTES 


Preliminary simulations were conducted to see the behavior of STARS and its 
output data. One of the test simulations crossed the equator to verify the system was not 
confined to the Northern Hemisphere. Another test simulation only used three UBs and 
LBs instead of four to see if the number of bounds could be altered. Once an 
understanding of how certain inputs generated invalid outputs was reached, systematic 
simulations were started. Three routes that are frequently traveled by USN ships were 
chosen: 

• Naval Station Norfolk, VA, to Naval Station Rota, Spain 

• Naval Station Pearl Harbor, HI, to Naval Station Yokosuka, Japan 

• Naval Station San Diego, CA, to Naval Station Guam 

All routes were run using environments for the period June 2010 to December 
2010. This period was chosen mostly due to the fact that environmental data were 
available starting in March 2010. Another reason for the selection of this timeframe was 
to have the opportunity to use environmental conditions over three of the four seasons: 
late summer, fall and winter. For each origin-destination pair two simulations were 
conducted each month, two weeks apart, with one simulation in the month of December 
for a total of thirteen simulations per origin-destination pair. The effect of an arrival date 
constraint was not explored in these experiments. The arrival date chosen for each 
origin-destination pair was the last day of the year, December 31, 2010. This arrival date 
was late enough that no simulation approached it; therefore, the arrival date did not 
constrain the routes. 

Tables 4 to 6 show the origin-destination pairs and the departure dates and arrival 
dates of each simulation. Figures 3 to 7 show Google Earth images that depict the 
different routes that were generated for selected departure dates. STARS produces routes 
for each environment unless it cannot find a feasible route that does not exceed the SOE 
limitations in the input file within that environment. If STARS cannot find a route that 
does not exceed the SOE limitations, it does not produce a route. Given there are sixteen 
environments, STARS can produce between zero and sixteen different routes for a single 

origin-destination pair for a specified departure date. In some cases, two or more 
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environments will produce identical routes. Figures 3-7 show Google Earth images that 
include all of the routes generated by STARS for the entire set of origin-destination pairs 
for each route. Figures 6 and 7 show the differences in the number of routes generated by 
STARS for two specified departure dates. Given there are sixteen environments plus 
“analysis,” the route shown in Figure 6 only generated one route from the environments 
plus the “analysis” route, yielding a total of two routes. In Figure 7, for each 
environment plus the “analysis,” a different route was generated, yielding a total of 
seventeen routes. 


Departure Date 

Arrival Date 

(MMDDYYYY) 

(MMDDYYYY) 

06012010 

06082010 

06152010 

06222010 

07012010 

07082010 

07152010 

07222010 

08012010 

08082010 

08152010 

08222010 

09012010 

09082010 

09152010 

09222010 

10012010 

10082010 

10152010 

10222010 

11012010 

11082010 

11152010 

11222010 

12012010 

12082010 


Table 4. Departure and Arrival Dates for the Routes Between Norfolk, 

VA and Rota, Spain 
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Departure Date 

Arrival Date 

(MMDDYYYY) 

(MMDDYYYY) 

06022010 

06102010 

06162010 

06242010 

07022010 

07102010 

07162010 

07242010 

08022010 

08102010 

08162010 

08242010 

09022010 

09102010 

09162010 

09242010 

10022010 

10102010 

10162010 

10242010 

11022010 

11102010 

11162010 

11242010 

12022010 

12102010 


Table 5. Departure and Arrival Dates for the Routes Between Pearl Harbor, HI and 

Yokosuka, Japan 


Departure Date 

Arrival Date 

(MMDDYYYY) 

(MMDDYYYY) 

06032010 

06142010 

06172010 

06282010 

07032010 

07142010 

07172010 

07282010 

08032010 

08172010 

08172010 

08282010 

09032010 

09172010 

09172010 

09282010 

10032010 

10142010 

10172010 

10282010 

11032010 

11142010 

11172010 

11282010 

12032010 

12142010 


Table 6. Departure and Arrival Dates for the Routes Between San Diego, CA and 

Naval Station Guam 
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Figure 3. Norfolk, VA to Rota, Spain (06012010), Showing STARS-generated 
Routes for Each Environment (From Google Earth, 2011) 



Figure 4. Pearl Harbor, HI, to Yokosuka, Japan (06022010), Showing STARS- 
generated Routes for Each Environment (From Google Earth, 2011) 
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Figure 5. San Diego, CA, to Naval Station Guam (06032010), Showing STARS- 
generated Routes for Each Environment (From Google Earth, 2011) 



Figure 6. Pearl Harbor, HI, to Yokosuka, Japan (08162010), Shows how 

Environments can be so Similar that the Same Route is Generated for 
Nearly all of the Environments (From Google Earth, 2011) 
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Figure 7. Seventeen Different Routes Generated by STARS for the Origin-destination 
Pair of Norfolk, VA, to Rota, Spain (12012010) (From Google Earth, 2011) 


17 




THIS PAGE INTENTIONALLY LEFT BLANK 


18 



III. ANALYSIS AND RESULTS 


A. OVERVIEW 

Every departure date for each origin-destination pair was analyzed to determine 
the perfonnance, the sensitivity and the robustness of the routes generated. The 
performance of the routes is measured as the amount of energy consumed and the 
distance traveled. The sensitivity of the routes reflects how far apart optimized routes 
are, measured in tenns of the distances among routes optimized with respect to different 
environmental forecasts. Robustness reflects how well the optimized route for a given 
origin-destination pair and start date stands up against the different environmental 
conditions, measured with respect to the safety of each route when environmental 
conditions are different from those used by the optimizer. 

B. PERFORMANCE ANALYSIS AND RESULTS 

The perfonnance of the routes is measured as the amount of energy consumed and 
the distance traveled as an average of the indicated number of optimized routes generated 
for each specified departure date, as shown in Tables 7-9. As discussed earlier, STARS 
did not generate an optimized route for every environment. For each origin-destination 
pair and departure date, M' is the number of ensemble members for which an optimized 
route was generated. These perfonnance tables show energy consumption (in kilo- 
horsepower-hours, khph) and distance (in nautical miles, nm) along the optimized route 
with average and standard deviation over the M' ensemble members. The average time 
traveled was approximately the same for each specified origin-destination pair plus or 
minus an hour, regardless of the departure date. The Norfolk to Rota route took 
179 hours. The Pearl Harbor to Yokosuka route took 189 hours, and the San Diego to 
Guam route took 270 hours. 
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For the Norfolk-Rota results, the average energy consumption and distances were 
similar during the earlier months. During the latter months, the energy consumption and 
distances increased, with the exception of 11012010. This reflects the generally more 
adverse environmental conditions during the later period, with the exception of 
11012010, reflecting variability even within the fall season. 


Origin- 

Destination 

Departure Date 
(MMDDYYYY) 

Total Energy 
Consumption 
(khph) 

Average (StDev) 

Total Distance 
(nm) Average 
(StDev) 

Number of 
Optimized 
Routes 
Generated 
(M’> 


06012010 

1853.3(10.2) 

3719.2(5.2) 

16 


06152010 

1864.7(12.8) 

3728.7(4.0) 

16 


07012010 

1823.0(9.0) 

3720.8(3.4) 

16 


07152010 

1848.8(8.5) 

3714.6(5.0) 

16 

Gg 

08012010 

1810.2(11.6) 

3714.3(1.9) 

16 

o 

(V 

08152010 

1867.3(9.1) 

3721.9(1.9) 

16 

HH 

09012010 

1826.5(7.0) 

3720.8(2.7) 

16 

a 

09152010 

1845.3(7.8) 

3718.2(5.7) 

16 

u 

© 

10012010 

1944.0(19.1) 

3735.0(1.6) 

16 

z 

10152010 

1932.6(33.9) 

3735.4(9.1) 

16 


11012010 

1890.9(29.2) 

3722.6(4.8) 

16 


11152010 

1975.2(44.8) 

3732.5(13.5) 

16 


12012010 

2197.1(93.9) 

3747.6(27.2) 

16 


12152010 

2016.6(42.6) 

3736.0(12.0) 

14 


Table 7. For Each Departure Date for the Norfolk-Rota Route, Total Horsepower 
and Distance Traveled on Optimized Route, Averaged Over all Ensemble 
Members That Produced an Optimized Route 


For the Pearl Harbor-Yokosuka results, the average energy consumption is similar 
for most of the departure dates. The largest energy requirements were in November and 
December. The average total distance traveled showed small differences of no more than 
50 mn over all departure dates. 
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Origin- 

Destination 

Departure Date 
(MMDDYYYY) 

Total Energy 
Consumption 
(khph) 

Average (StDev) 

Total Distance 
(nm) 

Average (StDev) 

Number of 
Optimized 
Routes 
Generated 
(M’) 


06022010 

2219.4(14.8) 

4079.1(0.0) 

16 


06162010 

2320.0(37.5) 

4091.2(0.1) 

15 


07022010 

2225.4(19.0) 

4106.3(10.5) 

15 

S3 

j3 

07162010 

2194.6(23.1) 

4087.5(5.2) 

16 

C/3 

O 

k> 

08022010 

2215.8(16.9) 

4083.9(3.7) 

15 

p-£ 

o 

> 

08162010 

2216.1(16.3) 

4068.7(0.5) 

16 


09022010 

2196.7(7.8) 

4103.8(2.6) 

16 

o 

pO 

09162010 

2250.1(9.2) 

4102.8(6.5) 

16 

Lm 

CS 

ffl 

10022010 

2232.1(52.5) 

4072.0(6.8) 

16 

Tn 

10162010 

2258.2(44.4) 

4071.0(6.7) 

16 

QJ 

Cu 

11022010 

2603.7(68.0) 

4082.5(7.9) 

15 


11162010 

2522.4(66.5) 

4082.9(14.1) 

16 


12022010 

2413.5(91.6) 

4078.4(10.4) 

16 


12162010 

2360.8(62.0) 

4056.9(7.2) 

16 


Table 8. For Each Departure Date for the Pearl Harbor-Yokosuka Route, Total 

Energy Consumption and Distance Traveled on Optimized Route, Averaged 
Over all Ensemble Members That Produced an Optimized Route 


For the San Diego-Guam results, note that there are no results shown from 
departure date 12172010, presumably because adverse conditions prevented STARS from 
finding an optimal route that satisfied constraints. As with the other origin-destination 
pairs, the average energy consumption was greater during the winter months. The 
highest average energy consumption was experienced during December. The average 
total distances traveled varied by +/- 20 mn. 
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Origin- 

Destination 

Departure Date 

(MMDDYYYY) 

Total Energy 

Consumption 

(khph) 

Average (StDev) 

Total Distance 

(nm) 

Average (StDev) 

Number of 

Optimized 

Routes 

Generated 

(M’) 


06032010 

3670.4(66.7) 

6065.4(0.0) 

14 


06172010 

3783.2(61.1) 

6054.7(0.0) 

15 


07032010 

3772.3(50.1) 

6059.7(0.0) 

16 


07172010 

3564.9(4.9) 

6053.1(0.0) 

14 

S 

08032010 

3560.5(14.7) 

6056.5(6.4) 

16 

s 

0 

08172010 

3523.2(9.5) 

6047.4(0.0) 

15 

o 

09032010 

3652.3(22.4) 

6068.4(3.4) 

15 

OX) 

QJ 

09172010 

3710.7(38.4) 

6057.9(0.0) 

16 

5 

10032010 

3742.6(34.5) 

6081.7(14.5) 

16 

Vl 

10172010 

3809.8(26.8) 

6079.7(6.1) 

16 


11032010 

4073.7(219.9) 

6044.4(0.0) 

15 


11172010 

3996.1(92.1) 

6062.5(0.0) 

16 


12032010 

4126.0(88.8) 

6047.8(0.0) 

16 


12172010 

NaN 

NaN 

NaN 


Table 9. For Each Departure Date for the San Diego-Guam Route, Total Energy 
Consumption and Distance Traveled on Optimized Route, Averaged Over 
all Ensemble Members That Produced an Optimized Route 


C. MEASURING SENSITIVITY 

The sensitivity metrics reflect how much the routes generated from a given origin- 
destination pair and start date differ from each other. Tables 10-18 summarize how far 
apart the routes optimized with respect to the 16 ensemble members lie. Specifically, 
tables give the average (over M' ensemble members) of the average and the maximum 
(over the route) great-circle (GC) distance between each ensemble member’s route and 
reference routes. The reference routes are: 

• The average of the optimized routes from the 16 ensemble members 
(Tables 10, 13 and 16); 

• The route optimized with respect to the analysis ( 11, 14, and 17); and 
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• The route optimized with respect to the average of the M —16 ensemble- 
member environments (Tables 12, 15, and 18). 

The distance between routes is measured in Tables 10-19 with both the average 

(Equation 1) and maximum over N — 2 waypoints, excluding the origin and destination 

of the GC distance between routes. The Routex environment used to optimize (£) is 

indicated in the row header, and j is the environment used to optimize reference route 

given in the column header. 

——r Yj [s reat - circle (*» (0» y\ O '),K U), y n (./'))) [ 1 ] 

^ n-2 


For Tables 10, 13, and 16, the results shown are averages (formula shown in Equation 2) 
and standard deviations over the number of ensemble members (M') for which STARS 
could generate an optimal route. 


1 1 

M'^Xn-2 



[ 2 ] 


Each table also shows the maximum GC distance between routes over the N waypoints, 
as in Equation 3 and, for Tables 10, 13, and 16, the average and standard deviation over 
all M' ensemble members are shown. 

max (great _ circle ( x* (i ), v* (i ), x* (j), v* (./))) [ 3 ] 

Each origin-destination pair has departure dates with noteworthy results. Tables 
10-12 show results for the Norfolk to Rota routes. Although the Norfolk to Rota routes 
show smallest maximal distances between routes, optimized routes for the same departure 
date may be more than 500 mn apart. Tables 10-12 show three departure dates for the 
Norfolk-Rota route that include an asterisk (*). Optimized routes for these departure 
dates, to include the analysis and the environment average routes, are shown in Figures 
8 - 10 . 
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There are four departure dates (07012010, 07152010, 11012010, and 12012010) 
for which the analysis-optimized route precisely coincides with the route optimized for 
the averaged-environments. Therefore, results for these two routes are the same and the 
distances between these two routes are zero. 




Route-Average GC Distance (nm) 
from 

Maximum GC Distance (nm) from 


Departure Date 
(MMDDYYYY) 

Analysis- 

optimized 

Route 

Averaged 

Environments- 

optimized 

Route 

Average of 
Member- 
optimized 
Route 

Analysis- 

optimized 

Route 

Averaged 

Environments- 

optimized 

Route 

Average of 
Member- 
optimized 
Route 

All Ensemble Members Average (StDev) 

06012010 

50.8(18.2) 

38.9(28.0) 

34.5(11.2) 

195.8(30.6) 

227.3(69.1) 

153.5(27.6) 

06152010 

108.6(7.0) 

12.4(17.0) 

15.4(12.4) 

332.6(33.0) 

184.3(55.2) 

165.1(42.7) 

07012010 

41.9(16.2) 

41.9(16.2) 

33.9(7.7) 

205.2(62.2) 

205.2(62.2) 

131.7(12.1) 

07152010 

13.9(25.6) 

13.9(25.6) 

20.6(14.8) 

121.9(45.8) 

121.9(45.8) 

99.8(26.5) 

08012010* 

30.6(15.1) 

8.4(13.9) 

12.4(9.3) 

239.9(75.2) 

194.5(72.9) 

158.2(47.6) 

08152010 

24.2(19.6) 

35.5(17.2) 

23.8(11.3) 

166.6(60.2) 

189.2(46.3) 

130.3(27.8) 

09012010 

47.7(12.5) 

28.4(27.6) 

33.2(12.9) 

227.2(44.0) 

389.8(148.2) 

290.7(78.5) 

09152010 

49.7(19.2) 

38.3(28.8) 

34.0(21.9) 

205.1(53.4) 

270.2(69.6) 

196.4(44.2) 

10012010* 

26.7(21.0) 

10.0(22.0) 

14.9(17.9) 

335.6(71.5) 

335.6(82.3) 

317.3(72.7) 

10152010 

72.8(27.1) 

59.4(20.5) 

57.5(11.4) 

360.2(74.4) 

295.9(52.0) 

243.0(39.4) 

11012010 

41.7(44.4) 

41.7(44.4) 

47.3(35.0) 

366.8(102.0) 

366.8(102.0) 

336.5(79.4) 

11152010 

89.8(61.3) 

107.2(49.8) 

76.4(22.5) 

463.7(152.0) 

375.4(99.8) 

295.8(53.8) 

12012010* 

120.1(79.2) 

120.1(79.2) 

104.5(24.8) 

503.9(187.1) 

503.9(187.1) 

380.6(79.8) 

12152010 

85.3(57.9) 

95.6(27.1) 

75.1(15.8) 

403.1(160.4) 

316.0(79.8) 

221.3(33.4) 


Table 10. For Each Departure Date for the Norfolk-Rota Route, Route Average GC 
Distance and Maximum GC Distance From Reference Routes (Summary 
Statistics Over Routes Optimized for Each Ensemble Member) 
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Route-Average GC Distance 
(nm) from 

Maximum GC Distance (nm) 
from 


Departure Date 
(MMDDYYYY) 

Analysis-optimized 

Route 

Averaged 
Environments- 
optimized Route 

Average of 
Member-optimized 
Route 

Analysis-optimized 

Route 

Averaged 
Environments- 
optimized Route 

Average of 
Member-optimized 
Route 


06012010 

0 

33.3 

42.8 

0 

159.9 

103.1 


06152010 

0 

109.6 

107.2 

0 

310.7 

297.1 


07012010 

0 

0 

34.5 

0 

0 

105.8 

W 

= 

07152010 

0 

0 

13.5 

0 

0 

24.4 

o 

P4 

08012010* 

0 

35.7 

30.4 

0 

194.5 

158.2 

a j 

N 

08152010 

0 

32.3 

17.5 

0 

77.2 

65.0 

1 

09012010 

0 

41.2 

34.5 

0 

227.2 

128.7 


09152010 

0 

33.3 

42.8 

0 

159.9 

103.1 

w 

10012010* 

0 

17.8 

25.0 

0 

49.8 

52.1 


10152010 

0 

73.2 

54.4 

0 

235.5 

154.1 

a 

c 

11012010 

0 

0 

18.7 

0 

0 

52.4 


11152010 

0 

157.2 

80.4 

0 

359.3 

186.3 


12012010* 

0 

0 

93.7 

0 

0 

175.2 


12152010 

0 

110.7 

67.4 

0 

316.0 

205.9 


Table 11. For Each Departure Date for the Norfolk-Rota Route, Route Average GC 
Distance and Maximum GC Distance From the Analysis-optimized Route 
and the Analysis-optimized Route, the Averaged Environments-optimized 
Route and the Average of Member-optimized Route 
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Averaged Environments- Optimized Route 


Route-Average GC Maximum GC Distance (nm) 

Distance (nm) from_from_ 



Table 12. For Each Departure Date for the Norfolk-Rota Route, Route Average GC 
Distance and Maximum GC Distance From the Averaged Environments- 
optimized Route and the Analysis-optimized Route, the Averaged 
Environments-optimized Route and the Average of Member-optimized 

Route 
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Figure 8. Optimized Routes Generated for Norfolk-Rota Route, 08012010 (From 

Google Earth, 2011) 
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Figure 9. Optimized Routes Generated for Norfolk-Rota Route, 10012010 (From 

Google Earth, 2011) 
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Figure 10. Optimized Routes Generated for Norfolk-Rota Route, 12012010 (From 

Google Earth, 2011) 

Tables 13-15 show sensitivity of routes from Pearl Flarbor to Yokosuka. This 
origin-destination pair shows the greatest sensitivity, with some optimized routes over 
5000 nm apart, for the same departure date. There are five departure dates (06162010, 
10022010, 10162010, 11022010 and 12162010) in which the analysis route and the 
ensemble-average-optimized route are the same. This does not necessarily mean that all 
routes are close to each other; rather, there are some ensemble members whose routes 
may be far (almost 780 nm away for 10162010) from the reference route. Different 
ensemble members may produce dramatically different routes. 
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All Ensemble Members Average(StDev) 



Table 13. For Each Departure Date for the Pearl Harbor-Yokosuka Route, Route 
Average GC Distance and Maximum GC Distance From the Average of all 
Ensemble Members and the Analysis-optimized Route, the Averaged 
Environments-optimized Route and the Average of Member-optimized 

Route 
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Analysis-Optimized Route 


Route-Average GC 


Maximum GC Distance 
from 


Departure Date 
(MMDDYYYY) 

Analysis- 
optimized Route 

Average 
Environments- 
optimized Route 

Average of 
Member- 
optimized Route 

Analysis- 
optimized Route 

Average 
Environments- 
optimized Route 

Average of 
Member- 
optimized Route 

06022010 

0 

NaN 

0 

0 

NaN 

0 

06162010 

0 

0 

122.3 

0 

0 

1134.4 

07022010* 

0 

49.4 

345.8 

0 

162.4 

6034.6 

07162010* 

0 

NaN 

207.3 

0 

NaN 

1261.2 

08022010 

0 

63.8 

61.1 

0 

216.4 

201.0 


08162010 

09022010 

09162010 

10022010 

10162010 

11022010 

11162010 

12022010 

12162010* 


17.2 
1.0 

53.3 

_ 0 _ 

_ 0 _ 

_ 0 _ 

NaN 

257.5 

0 


17.1 

88.7 

227.3 

129.2 
435.6 

274.8 

132.2 

133.9 
3.8 


94.3 

31.4 

144.5 

_0_ 

_0_ 

_0_ 

NaN 

647.5 
0 


90.5 

2372.1 

5725.3 

1208.5 

3113.3 

2234.6 
1068.0 

1144.3 
18.7 


Table 14. For Each Departure Date for the Pearl Harbor-Yokosuka Route, Route 
Average GC Distance and Maximum GC Distance From the Analysis- 
optimized Route and the Analysis-optimized Route, the Averaged 
Environments-optimized Route and the Average of Member-optimized 

Route 
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Average Environments-Optimized Route 


Route-Average GC Distance Maximum GC Distance (nm) 
_(nm) from_from_ 



Table 15. For Each Departure Date for the Pearl Harbor-Yokosuka Route, Route 
Average GC Distance and Maximum GC Distance From the Averaged 
Environments-optimized Route and the Analysis-optimized Route, the 
Averaged Environments-optimized Route and the Average of Member- 

optimized Route 
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The fields in Tables 14 and 15 with “NaN” as the cell value mean that the 
optimized route used for comparison was not generated. Tables 13-15 show three 
departure dates for the Pearl Harbor-Yokosuka route that include an asterisk. Optimized 
routes for these departure dates, to include the analysis and the environment average 
routes, are shown in Figures 11-13. 



Figure 11. Optimized Routes Generated for Pearl Harbor-Yokosuka Route, 06162010 

(From Google Earth, 2011) 
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Figure 12. Optimized Routes Generated for Pearl Harbor-Yokosuka Route, 07162010 

(From Google Earth, 2011) 
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Figure 13. Optimized Routes Generated for Pearl Harbor-Yokosuka Route, 12162010 

(From Google Earth, 2011) 
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Tables 16-18 show sensitivity results for routes from San Diego to Guam. There 
are only four departure dates (08032010, 09032010, 10032010, and 10172010) with any 
sensitivity. However, on these four dates, the optimized routes are very far apart—almost 
800 nm for 10032010 departure date. 




Route-Average GC Distance (nm) 
from 

Maximum GC Distance from 


Departure Date 
(MMDDYYYY) 

Analysis-optimized 

Route 

Average 
Environments- 
optimized Route 

Average of Member- 
optimized Route 

Analysis-optimized 

Route 

Average 
Environments- 
optimized Route 

Average of Member- 
optimized Route 

All Ensemble Members Average (StDev) 

06032010 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

06172010 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

07032010 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

07172010 

NaN(NaN) 

0.0(0.0) 

0.0(0.0) 

NaN(NaN) 

0.0(0.0) 

0.0(0.0) 

08032010* 

71.3(31.6) 

69.2(31.8) 

200.4(10.3) 

336.0(82.6) 

336.0(85.4) 

6174.8(26.6) 

08172010 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

09032010* 

18.9(73.2) 

18.9(73.2) 

35.3(63.5) 

756.2(195.2) 

756.2(195.2) 

705.9(169.3) 

09172010 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

10032010* 

203.6(67.9) 

133.2(96.2) 

206.7(36.7) 

604.6(111.6) 

790.2(262.5) 

2259.9(77.4) 

10172010* 

195.0(20.1) 

47.8(80.9) 

176.6(48.3) 

512.8(86.2) 

555.9(206.2) 

5029.9(31.6) 

11032010 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

11172010 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

12032010 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

0.0(0.0) 

12172010 

NaN(NaN) 

NaN(NaN) 

NaN(NaN) 

NaN(NaN) 

NaN(NaN) 

NaN(NaN) 


Table 16. For Each Departure Date for the San Diego-Guam Route, Route Average 
GC Distance and Maximum GC Distance From the Average of all 
Ensemble Members and the Analysis-optimized Route, the Averaged 
Environments-optimized Route and the Average of Member-optimized 

Route 
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Optimized Route 


Route-Average GC Distance 

_(nm) from_Maximum GC Distance from 



Table 17. For Each Departure Date for the San Diego-Guam Route, Route Average 
GC Distance and Maximum GC Distance From the Average of the 
Analysis-optimized Route and the Analysis-optimized Route, the Averaged 
Environments-optimized Route and the Average of Member-optimized 

Route 
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Average Environments-Optimized Route 


Route-Average GC Distance 

_(nm) from_Maximum GC Distance from 



Table 18. For Each Departure Date for the San Diego-Guam Route, Route Average 
GC Distance and Maximum GC Distance From the Average Environments- 
optimized Route and the Analysis-optimized Route, the Averaged 
Environments-optimized Route and the Average of Member-optimized 

Route 
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Figure 14. Optimized Routes Generated for San Diego-Guam Route, 08032010 (From 

Google Earth, 2011) 
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Figure 15. Optimized Routes Generated for San Diego-Guam Route, 09032010 (From 

Google Earth, 2011) 
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Figure 16. Optimized Routes Generated for San Diego-Guam Route, 10032010 (From 

Google Earth, 2011) 
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Figure 17. Optimized Routes Generated for San Diego-Guam Route, 10032010 

(From Google Earth, 2011) 


In some cases, the average of the all ensemble member routes produced larger 
average GC distances from the member-optimized routes than the analysis-optimized 
route or the environments-optimized routes. This is counterintuitive. The reason is that 
the average route consists of the Cartesian average position of each waypoint, while the 
distances between routes are calculated using GC distances between waypoints. Overall, 
the Norfolk-Rota and Pearl Harbor-Yokosuka routes are more sensitive to the 
environmental conditions during the time period indicated. The San Diego-Guam routes 
are the least sensitive to the environmental conditions. With only four departure dates 
with non-zero results, this means that the environmental conditions on the other departure 
dates have no effect on the route generated. Regardless of the environmental conditions, 
the routes generated are optimal in any environmental condition. 
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D. MEASURING ROBUSTNESS 


The robustness metrics reflect how well the optimized route for a given 
environment origin-destination pair and start date measures up with respect to the safety 
thresholds when evaluated within other forecast environmental conditions. 

In addition to the variables generated by STARS, the following variables were 
calculated as they will be used in measuring route safety: 


relative _ wind ( w n , j ) 

Wind speed experienced by ship while in 
motion. 

seas(w n J) 

The maximum 

of sig _ wave _ height ( w n , /) 

sea _ height ( w n , j) 
and s weU height(w n ,j) 


Relative wind was calculated according to Equation 4. 

yj wind _ speed (w n ,jf + s n + 2 ( wind _ speed ( w n , y)H x cos C n [4] 

The safety thresholds utilized in STARS, as depicted in Appendix A, are 35 knots 
of wind speed, 65 knots of relative wind speed, and 25-foot seas (head, beam and 
following). The STARS output did not provide head, beam, and following seas. In 
evaluating the safety of routes, the operational threshold is 25-foot seas, 
using seas (w„,y). The notation condition ( w n , j ) will be used to generically represent 

any of the three METOC conditions used to measure safety, i.e., wind _speed ( w n , jj , 

relative_wind ( w n ,j ) or seas ( w n , j ). 

Tables 19-27 show summary statistics for safety-related METOC conditions and 
energy consumption relative to the same environment used for optimization. For each 
origin-destination pair and start date, these tables give the mean and maximum 
environmental conditions that would be experienced by the optimized route if the forecast 
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conditions hold, as well as the number of safety threshold exceedances experienced (as 
shown below in Equations 5-7, respectively), and the average amount of the exceedance, 
when it occurred. 


N 


N 

y, condition ( w* (/),/j 

n =1 


max condition[w* n (/j, /j 


[5] 

[ 6 ] 


1 N 

exceedances = — '^X n (i,j) 

N „=i 


x.(uh 


1 if condilioni vv'(;)./| > condition threshold 
0 otheiwise 


[7] 


Tables 19-21 show the overall safety of the optimized routes for each origin- 
destination pair, evaluated with respect to their own environment. The fonnulas for the 
mean and maximum values are shown in Equations 8 and 9, respectively. In the top third 
of Tables 19-21, mean and standard deviation with respect to the environments for which 
the routes were optimized ( i ) are shown. 


N 


^ condition ( vv* (/),/) 

n =1 


max 

n=l,...,N 


| condi tion (vv* (i),ijj 


[ 8 ] 

[9] 
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Absolute Winds, kts 

(35kt Threshold) 

Relative Winds, <ts(65kt Threshold) 

Seas, ft (25ft Thresh old) 


Departure 

Date 

(MMDDYYYY) 

Energy 

Consumption 

(khph) 

Mean 

X 

RE 

2 

Exceedances 

(Number) 

Exceedances 

(kts) 

C 

ra 

a 

2 

X 

IS 

s 

Exceedances 

(Number) 

I/I 

8 

c 

rs 

■8 

S a 
£ 2. 

Mean 

X 

ra 

s 

Exceedances 

(Number) 

Exceedances 

(ft) 


06012010 

1853.3(10.2) 

16.4(1.64) 

26.5(1.77) 

0.0(0.00) 

NaN 

346(1.24) 

53.7(2.34) 

0.0(0.00) 

NaN 

7.7(0.50) 

11.1(0.74) 

0.0(0.00) 

NaN 


06152010 

18617(12.8) 

10.6(1.50) 

240(3.02) 

0.0(0.00) 

NaN 

26.1(168) 

39.2(2.75) 

0.0(0.00) 

NaN 

5.4(0.60) 

10.5(0.90) 

0.0(0.00) 

NaN 

a 

07012010 

1823.0(9.0) 

12.9(1.07) 

27.1(3,54) 

0.0(0.00) 

NaN 

26.6(0.94) 

43.0(2.52) 

0.0(0.00) 

NaN 

47(0.35) 

10.8(108) 

0.0(0.00) 

NaN 


07152010 

1848.8(8.5) 

141(1.32) 

22.2(1.01) 

0.0(0.00) 

NaN 

33.6(1.26) 

43.3(1.96) 

0.0(0.00) 

NaN 

6.2(0.56) 

9.310.66) 

0.0(0.00) 

NaN 

c 

& 

08012010 

1810.2(11.6) 

12.2(1.65) 

26.8(3.88) 

0.0(0.00) 

NaN 

25.0(1.16) 

44.9(4.34) 

0.0(0.00) 

NaN 

3.7(0.42) 

9.11162) 

0.0(0.00) 

NaN 

2 

08152010 

1867.3(9.1) 

9.7(1.19) 

18.3(2.14) 

0.0(0.00) 

NaN 

26.4'137) 

41.4(2.67) 

0.0(0.00) 

NaN 

5.2(0.33) 

7.910.62) 

0.0(0.00) 

NaN 

c 

© 

09012010 

1826.5(7.0) 

9.9(0.78) 

19.6(2.12) 

0.0(0.00) 

NaN 

28.4125) 

45.3(2.43) 

0.0(0.00) 

NaN 

49(0.33) 

13.9(131) 

0.0(0.00) 

NaN 

£ 

09152010 

1845.3(7.8) 

11.4(2.02) 

28.0(3.43) 

0.0(0.00) 

NaN 

29.0(2.15) 

53.6(6.38) 

0.0(0.00) 

NaN 

5.8(0.62) 

10.21106) 

0.0(0.00) 

NaN 

s 

10012010 

19440(19.1) 

15.6(1.11) 

27.2(2.06) 

0.0(0.00) 

NaN 

32.0(1.64) 

48.9(2.36) 

0.0(0.00) 

NaN 

9.1(0.68! 

17.8(154) 

0.0(0.00) 

NaN 

-t 

-□ 

10152010 

1932.6(33.9) 

18.9(1.22) 

348(2.34) 

0.0(0.00) 

NaN 

33.9(1.63) 

59.7(3.43) 

0.0(0.00) 

NaN 

10.9(0.83) 

19.5(2.22) 

0.0(0.00) 

NaN 

E 

as 

11012010 

1890.9(29.2) 

16.4(1.20) 

35.9(444) 

0.1(0.25) 

0.2(0.41) 

27.2(1.35) 

46.9(2.74) 

0.0(0.00) 

NaN 

7.5(0.53) 

15.112.06) 

0.0(0.00) 

NaN 

c 

LU 

11152010 

1975.2(44.8) 

13.3(1.95) 

36.7(477) 

0.2(0.40) 

1.0(0.70) 

24.1(1.78) 

51.7(6.47) 

0.0(0.00) 

NaN 

7.1(0.83) 

17.5(2.63) 

0.0(0.00) 

NaN 


12012010 

2197.1(93.9) 

18.3(1.95) 

31.6(1.49) 

0.0(0.00) 

NaN 

29.7(2.12) 

57.2(4.63) 

0.0(0.00) 

NaN 

11.7(1.30) 

247(3.79) 

0.0(0.00) 

NaN 


12152010 

2016.6(42.6) 

17.1(6.95) 

35.7(10.80) 

0.1(0.25) 

0.2(0.41) 

30.7(12.15) 

56.6(17.81) 

0.0(0.00) 

NaN 

10.0(3.99) 

20.8(6.48) 

0.0(0.00) 

NaN 


06012010 

1853.3(10.2) 

15.6 

21.3 

0 

NaN 

33.5 

43,0 

0 

NaN 

6.4 

8.7 

0 

NaN 


06152010 

18647(12.8) 

14.7 

22.1 

0 

NaN 

29.7 

40.0 

0 

NaN 

4.8 

8.2 

0 

NaN 


07012010 

1823.0(9.0) 

11.7 

18.4 

0 

NaN 

25.5 

36.5 

0 

NaN 

4.4 

6.8 

0 

NaN 

© 

07152010 

1848.8(8.5) 

142 

23.5 

0 

NaN 

33.2 

47.8 

0 

NaN 

6.2 

8.6 

0 

NaN 

g 

08012010 

1810.2(11.6) 

13.2 

23.7 

0 

NaN 

26.5 

42.1 

0 

NaN 

4.2 

8.1 

0 

NaN 

T3 

© 

08152010 

1867.3(9.1) 

10.8 

16.2 

0 

NaN 

28.3 

40.1 

0 

NaN 

5.4 

7.3 

0 

NaN 

£ 

09012010 

1826.5(7.0) 

12.2 

21.6 

0 

NaN 

31.0 

46.4 

0 

NaN 

6.1 

11.0 

0 

NaN 

a. 

09152010 

1845.3(7.8) 

12.9 

18.9 

0 

NaN 

29.9 

43.8 

0 

NaN 

6.5 

8.4 

0 

NaN 

.B 

10012010 

19440(19.1) 

140 

25.8 

0 

NaN 

28.8 

43.6 

0 

NaN 

9.3 

14.7 

0 

NaN 

j»> 

10152010 

1932.6(33.9) 

18.0 

27.7 

0 

NaN 

34.2 

52.6 

0 

NaN 

10.3 

14.3 

0 

NaN 

c 

«£ 

11012010 

1890.9(29.2) 

17.8 

243 

0 

NaN 

28.0 

42.1 

0 

NaN 

8.5 

11.5 

0 

NaN 


11152010 

1975.2(44.8) 

13.8 

30.0 

0 

NaN 

21.9 

33.4 

0 

NaN 

6.1 

ll.fi 

0 

NaN 


12012010 

2197.1(93.9) 

22.1 

30.7 

0 

NaN 

30.1 

46.3 

0 

NaN 

145 

25.9 

1 

0.85 


12152010 

2016.6(42.6) 

19.6 

31.5 

0 

NaN 

32.7 

53.1 

0 

NaN 

11.5 

18.0 

0 

NaN 


06012010 

1853.3(10,2) 

16.4 

21.5 

0 

NaN 

35.4 

46.3 

0 

NaN 

7.7 

9.9 

0 

NaN 

a 

'06152010 

18647(12.8) 

10.1 

15.4 

0 

NaN 

26.8 

32.2 

0 

NaN 

5.4 

8.5 

0 

NaN 

I 

r 07012010 

1823.0(9.0) 

12.1 

15.7 

0 

NaN 

26.3 

38.3 

0 

NaN 

4.5 

7.3 

0 

NaN 

■8 

'07152010 

1848.8(8.5) 

13.7 

19.2 

0 

NaN 

34.1 

46.6 

0 

NaN 

6.2 

7.6 

0 

NaN 

£ 

r 08012010 

1810.2(11.6) 

11.4 

16.3 

0 

NaN 

25.4 

40.8 

0 

NaN 

3.5 

5.7 

0 

NaN 

Q- 

r 08152010 

1867.3(9.1) 

9.4 

12.8 

0 

NaN 

28.0 

50.9 

0 

NaN 

5.3 

6.5 

0 

NaN 

M 

'09012010 

1826.5(7.0) 

9.7 

15.3 

0 

NaN 

28.8 

36.4 

0 

NaN 

5.0 

10.8 

Q 

NaN 

© 

E 

r 09152010 

1845.3(7.8) 

11.3 

15.9 

0 

NaN 

29.5 

40.3 

0 

NaN 

6.0 

7.4 

0 

NaN 

= 

£ 

r 10012010 

19440(19.1) 

145 

20.0 

0 

NaN 

33.1 

42.8 

0 

NaN 

9.0 

14.0 

0 

NaN 

> 

r 10152010 

1932.6(33.9) 

18.1 

31.5 

0 

NaN 

31.9 

50.2 

0 

NaN 

10.8 

15.4 

0 

NaN 

LU 

r 11012010 

1890.9(29.2) 

16.8 

23.2 

0 

NaN 

26.5 

41.0 

0 

NaN 

7.5 

9.9 

0 

NaN 

3C 

E 

'11152010 

1975.2(44.8) 

11.4 

210 

0 

NaN 

23.5 

38.6 

0 

NaN 

7.2 

10.9 

0 

NaN 

> 

r 12012010 

2197. K 93.9) 

19.0 

28.9 

0 

NaN 

29.7 

47.3 

0 

NaN 

13.2 

18.5 

0 

NaN 


'12152010 

2016.6(42.6) 

19.5 

27.2 

0 

NaN 

37.8 

50.5 

0 

NaN 

11.3 

16.2 

0 

NaN 


Table 19. Robustness of the Optimized Routes for Norfolk to Rota 
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Absolute Winds, kts (35kt Threshold) 

Relative Winds, kts (65kt Threshold) 

Seas, ft (25ft Threshold) 


Departure 

Date 

(MMDDYYYY) 

Energy 

Consumption 

(khph) 

Mean 

X 

TO 

Exceedances 

(Number) 

Exceedances 

(kts) 

Mean 

X 

TO 

Exceedances 

(Number) 

Exceedances 

(kts) 

Mean 

X 

ro 

Exceedances 

(Number) 

Exceedances 

(«) 

| All Ensemble Members Mean (StDev) | 

06022010 

2219.4(14.8) 

12.9(0.82) 

32.4(3.80) 

0.0(0.00) 

NaN 

28.0(1.19) 

49.4(2.69) 

0.0(0.00) 

NaN 

5.0(0.34) 

11.7(1.30) 

0.0(0.00) 

NaN 

06162010 

2320.0(37.5) 

14.4(4.05) 

30.0(6.80) 

0.0(0.00) 

NaN 

23.3(6.25) 

39.4(9.46) 

0.0(0.00) 

NaN 

6.4(1.80) 

14.1(3.09) 

0.0(0.00) 

NaN 

07022010 

2225.4(19.0) 

10.9(3.14) 

22.2(4.84) 

0.0(0.00) 

NaN 

24.2(6.54) 

42.7(10.27) 

0.0(0.00) 

NaN 

4.5(1.25) 

9.9(2.16) 

0.0(0.00) 

NaN 

07162010 

2194.6(23.1) 

11.1(4.78) 

26.6(7.54) 

0.0(0.00) 

NaN 

24.0(9.60) 

43.2(13.93) 

0.0(0.00) 

NaN 

3.7(1.59) 

10.6(2.70) 

0.0(0.00) 

NaN 

08022010 

2215.8(16.9) 

11.4(3.58) 

18.8(5.04) 

0.0(0.00) 

NaN 

26.4(7.54) 

45.9(11.75) 

0.0(0.00) 

NaN 

4.2(1.28) 

7.1(1.86) 

0.0(0.00) 

NaN 

08162010 

2216.1(16.3) 

12.3(0.81) 

24.5(2.16) 

0.0(0.00) 

NaN 

29.4(1.06) 

43.9(1.34) 

0.0(0.00) 

NaN 

4.8(0.57) 

10.9(1.52) 

0.0(0.00) 

NaN 

09022010 

2196.7(7.8) 

11.0(0.66) 

19.9(2.04) 

0.0(0.00) 

NaN 

28.4(0.90) 

44.3(1.90) 

0.0(0.00) 

NaN 

4.0(0.28) 

7.5(0.55) 

0.0(0.00) 

NaN 

09162010 

2250.1(9.2) 

9.8(0.70) 

18.2(0.62) 

0.0(0.00) 

NaN 

27.4(0.81) 

43.1(0.66) 

0.0(0.00) 

NaN 

5.1(0.20) 

9.1(0.59) 

0.0(0.00) 

NaN 

10022010 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

10162010 

2258.2(44.4) 

14.1(1.16) 

30.6(3.95) 

0.0(0.00) 

NaN 

29.3(0.96) 

51.6(4.36) 

0.0(0.00) 

NaN 

7.1(0.96) 

18.1(3.09) 

0.0(0.00) 

NaN 

11022010 

2603.7(68.0) 

15.2(4.16) 

33.6(7.66) 

0.0(0.00) 

NaN 

27.9(7.51) 

51.0(11.37) 

0.0(0.00) 

NaN 

10.2(2.80) 

17.4(3.98) 

0.0(0.00) 

NaN 

11162010 

2522.4(66.5) 

16.5(2.24) 

43.5(7.14) 

0.7(1.40) 

4.0(2.59) 

30.5(2.67) 

66.7(8.98) 

0.1(0.25) 

** 

12.0(1.20) 

27.7(3.76) 

0.1(0.50) 

** 

12022010 

2413.5(91.6) 

22.3(7.13) 

47.5(11.33) 

1.8(2.21) 

1.7(5.59) 

34.6(10.38) 

66.0(15.60) 

0.1(0.25) 

0.5(0.73) 

12.4(3.96) 

26.8(6.83) 

0.5(1.19) 

1.0(0.71) 

12162010 

2360.8(62.0) 

17.4(4.85) 

36.0(7.82) 

0.1(0.25) 

0.5(0.71) 

29.5(8.10) 

53.2(11.61) 

0.0(0.00) 

NaN 

9.4(2.58) 

15.3(3.61) 

0.0(0.00) 

NaN 

| Analysis-Optimized Route | 

06022010 

2219.4(14.8) 

14.1 

28.3 

0 

NaN 

28.1 

47.5 

0 

NaN 

5.8 

10.5 

0 

NaN 

06162010 

2320.0(37.5) 

16.8 

24.6 

0 

NaN 

24.5 

38.3 

0 

NaN 

7.0 

10.0 

0 

NaN 

07022010 

2225.4(19.0) 

12.2 

20.8 

0 

NaN 

26.1 

41.7 

0 

NaN 

4.5 

6.9 

0 

NaN 

07162010 

2194.6(23.1) 

NaN 

NaN 

NaN 

NaN 

23.5 

31.1 

0 

NaN 

NaN 

NaN 

NaN 

NaN 

08022010 

2215.8(16.9) 

11.7 

18.2 

0 

NaN 

28.9 

38.8 

0 

NaN 

5.1 

7.2 

0 

NaN 

08162010 

2216.1(16.3) 

13.9 

19.2 

0 

NaN 

30.1 

40.0 

0 

NaN 

5.2 

10.8 

0 

NaN 

09022010 

2196.7(7.8) 

10.6 

13.4 

0 

NaN 

28.8 

36.6 

0 

NaN 

4.3 

6.8 

0 

NaN 

09162010 

2250.1(9.2) 

10.3 

17.5 

0 

NaN 

28.7 

40.0 

0 

NaN 

5.8 

8.1 

0 

NaN 

10022010 

2232.1(52.5) 

15.3 

21.8 

0 

NaN 

27.0 

37.6 

0 

NaN 

6.7 

11.6 

0 

NaN 

10162010 

2258.2(44.4) 

NaN 

NaN 

NaN 

NaN 

23.7 

26.0 

0 

NaN 

NaN 

NaN 

NaN 

NaN 

11022010 

2603.7(68.0) 

16.6 

28.4 

0.0 

NaN 

25.6 

42.9 

0 

NaN 

10.7 

16.0 

0 

NaN 

11162010 

2522.4(66.5) 

NaN 

NaN 

NaN 

NaN 

23.8 

29.5 

0 

NaN 

NaN 

NaN 

NaN 

NaN 

12022010 

2413.5(91.6) 

24.1 

36.8 

3 

1.2 

35.6 

59.0 

0 

NaN 

13.6 

22.7 

0 

NaN 

12162010 

2360.8(62.0) 

18.7 

27.1 

0 

NaN 

33.2 

47.0 

0 

NaN 

9.9 

13.7 

0 

NaN 

| Average-Environment-Optimized Route | 

06022010 

2219.4(14.8) 

11.8 

15.2 

0 

NaN 

28.7 

40.4 

0 

NaN 

4.7 

6.8 

0 

NaN 

06162010 

2320.0(37.5) 

14.8 

23.3 

0 

NaN 

29.5 

43.0 

0 

NaN 

6.7 

10.0 

0 

NaN 

07022010 

2225.4(19.0) 

10.5 

17.2 

0 

NaN 

26.7 

41.5 

0 

NaN 

4.5 

7.5 

0 

NaN 

07162010 

2194.6(23.1) 

12.6 

15.1 

0 

NaN 

28.4 

41.3 

0 

NaN 

4.2 

4.8 

0 

NaN 

08022010 

2215.8(16.9) 

11.7 

17.3 

0 

NaN 

29.1 

44.6 

0 

NaN 

4.4 

6.5 

0 

NaN 

08162010 

2216.1(16.3) 

11.7 

18.8 

0 

NaN 

29.8 

38.7 

0 

NaN 

4.8 

8.8 

0 

NaN 

09022010 

2196.7(7.8) 

10.5 

14.2 

0 

NaN 

26.9 

36.7 

0 

NaN 

3.9 

5.7 

0 

NaN 

09162010 

2250.1(9.2) 

9.2 

17.1 

0 

NaN 

27.0 

42.1 

0 

NaN 

4.9 

7.4 

0 

NaN 

10022010 

2232.1(52.5) 

13.1 

18.8 

0 

NaN 

27.2 

38.3 

0 

NaN 

5.9 

8.3 

0 

NaN 

10162010 

2258.2(44.4) 

13.6 

18.6 

0 

NaN 

28.8 

42.8 

0 

NaN 

7.3 

13.2 

0 

NaN 

11022010 

2603.7(68.0) 

15.6 

22.6 

0 

NaN 

29.9 

42.2 

0 

NaN 

10.9 

15.8 

0 

NaN 

11162010 

2522.4(66.5) 

15.6 

26.7 

0 

NaN 

30.6 

41.2 

0 

NaN 

11.8 

15.1 

0 

NaN 

12022010 

2413.5(91.6) 

17.5 

22.7 

0 

NaN 

33.1 

45.0 

0 

NaN 

11.6 

17.6 

0 

NaN 

12162010 

2360.8(62.0) 

18.6 

27.1 

0 

NaN 

33.1 

43.7 

0 

NaN 

10.2 

14.3 

0 

NaN 


Table 20. Robustness of the Optimized Routes for Pearl Harbor to Yokosuka 
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Absolute Winds, Ids (35kt Threshold) 

Relative Winds, Ids (65ktThreshold) 

Seas, ft (25ftThreshold) 


Departure 

Date 

(MMDDYYYY) 

Mean 

Max 

Exceedances 

(Number) 

8 

C 

TJ 

<D 

8 a 

x 52 

LU ~ 

Mean 

Max 

Exceedances 

(Number) 

Exceedances 

(kts| 

Mean 

Max 

Exceedances 

(Number) 

Exceedances 

(ft) 

All Ensemble Members Mean (StDev) 

06032010 

11.1(4.49] 

32.7(9.35) 

0.0(0.00) 

NaN 

22.4(8.76) 

54.1(15.45) 

0.0(0.00) 

NaN 

5.9(249) 

21.3(4.97) 

0.0(0.00) 

NaN 

06172010 

17.1(4.76] 

31.0(6.94) 

0.0(0.00) 

NaN 

21.2(5.68) 

49.0(10.78) 

0.0(0.00) 

NaN 

8.1(227) 

14.5(3.37) 

0.0(0.00) 

NaN 

07032010 

17.6(1.26] 

34.2(3.00) 

0.0(0.00) 

NaN 

26.3(0.57) 

47.8(0.00) 

0.0(0.00) 

NaN 

8.6(0.69) 

15.8(1.24) 

0.0(000) 

NaN 

07172010 

9.5(3.81] 

21.2(6.76) 

0.0(0.00) 

NaN 

22.1(8.66) 

39.7(12.51) 

0.0(0.00) 

NaN 

3.8(151) 

7.9(2.30) 

0.0(000) 

NaN 

08032010 

11.6(1.01] 

18.3(1.32) 

0.0(0.00) 

NaN 

23.5(0.76) 

37.3(226) 

0.0(0.00) 

NaN 

4.8(0.29) 

7.7(0.72) 

0.0(000) 

NaN 

08172010 

8.0(230] 

16.6(3.82) 

0.0(0.00) 

NaN 

23.5(6.31) 

39.5(9.43) 

0.0(0.00) 

NaN 

3.4(0.93) 

6.5(1.49) 

0.0(000) 

NaN 

09032010 

13.0(3.71) 

23.7(5.55) 

0.0(0.00) 

NaN 

222(6.01) 

49.6(1264) 

0.0(0.00) 

NaN 

6.0(168) 

11.8(2.46) 

0.0(000) 

NaN 

09172010 

13.5(1.21] 

24.9(1.85) 

0.0(0.00) 

NaN 

23.5(0.79) 

42.5(250) 

0.0(0.00) 

NaN 

7.5(0.76) 

14.7(1.03) 

0.0(000) 

NaN 

10032010 

15.1(1.39] 

34.4(3.57) 

0.0(0.00) 

NaN 

26.1(0.97) 

56.6(5.93) 

0.0(0.00) 

NaN 

8.3(0.55) 

16.0(1.83) 

0.0(000) 

NaN 

10172010 

15.3(1.54] 

30.6(3.05) 

0.0(0.00) 

NaN 

26.1(0.85) 

52.0(4.19) 

0.0(0.00) 

NaN 

9.0(0.47) 

16.5(1.11) 

0.0(000) 

NaN 

11032010 

14.3(5.99] 

38.2(11.91) 

0.2(0.41) 

1.K1.31) 

21.2(8.31) 

43.6(13.47) 

0.0(0.00) 

NaN 

11.1(4.60) 

23.8(6.71) 

0.0(000) 

NaN 

11172010 

16.6(0.56] 

37.9(3.08) 

0.3(058) 

1.4(1.33) 

26.0(0.91) 

52.7(414) 

0.0(0.00) 

NaN 

11.9(0.57) 

20.8(1.11) 

0.0(000) 

NaN 

12032010 

21.7(1.81] 

56.9(6.43) 

1.9(182) 

41(2.72) 

30.6(1.51) 

63.8(3.91) 

0.0(0.00) 

NaN 

14.4(1.04) 

31.5(3.01) 

2.6(1.79) 

0.4(174) 

12172010 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

Analysis-Optimized Route 

06032010 

12.8 

20.3 

0 

NaN 

26.8 

41.1 

0 

NaN 

5.8 

10.8 

Q 

NaN 

06172010 

16.7 

39.1 

1 

4.1 

23.0 

51.2 

0 

NaN 

7.6 

14.6 

0 

NaN 

07032010 

19.3 

32.7 

0 

NaN 

27.3 

47.8 

0 

NaN 

8.4 

13.8 

0 

NaN 

07172010 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

08032010 

11.7 

15.5 

0 

NaN 

25.0 

37.9 

0 

NaN 

4.6 

6.5 

0 

NaN 

08172010 

NaN 

NaN 

NaN 

NaN 

24.7 

29.0 

0 

NaN 

NaN 

NaN 

NaN 

NaN 

09032010 

13.3 

24.1 

0 

NaN 

23.7 

49.6 

0 

NaN 

5.6 

9.4 

0 

NaN 

09172010 

14.6 

26.2 

0 

NaN 

24.8 

36.3 

0 

NaN 

7.9 

11.9 

0 

NaN 

10032010 

13.4 

27.4 

0 

NaN 

26.7 

42.0 

0 

NaN 

8.0 

11.8 

a 

NaN 

10172010 

15.9 

30.6 

0 

NaN 

26.0 

38.7 

0 

NaN 

8.5 

12.7 

0 

NaN 

11032010 

19.3 

28.2 

0 

NaN 

25.2 

48.2 

0 

NaN 

12.2 

16.2 

0 

NaN 

11172010 

20.5 

40.5 

2 

5.5 

29.5 

57.4 

0 

NaN 

15.3 

31.5 

4 

1.4 

12032010 

16.0 

37.1 

1 

2.1 

26.3 

54.1 

0 

NaN 

11.1 

18.8 

0 

NaN 

12172010 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

5 

3 

O 

cc 

o> 

1 

a. 

a 

6 
E 

C, 

o 

"> 

c. 

UJ 

i 

cno 

£ 

< 

06032010 

11.6 

18.9 

0 

NaN 

25.3 

37.4 

0 

NaN 

6.6 

9.9 

0 

NaN 

06172010 

17.2 

243 

0 

NaN 

28.4 

41.9 

0 

NaN 

8.3 

11.1 

0 

NaN 

07032010 

16.6 

23.7 

0 

NaN 

26.2 

47.8 

0 

NaN 

8.3 

12.7 

0 

NaN 

07172010 

10.5 

19.5 

0 

NaN 

25.5 

36.0 

0 

NaN 

4.4 

6.2 

0 

NaN 

08032010 

11.8 

16.0 

0 

NaN 

22.5 

34.3 

0 

NaN 

4.9 

6.5 

0 

NaN 

08172010 

8.2 

13.5 

0 

NaN 

25.4 

36.1 

0 

NaN 

3.6 

5.3 

0 

NaN 

09032010 

13.4 

20.3 

0 

NaN 

23.6 

49.6 

0 

NaN 

6.4 

8.4 

0 

NaN 

09172010 

12.3 

17.5 

0 

NaN 

22.3 

341 

0 

NaN 

7.2 

12.7 

0 

NaN 

10032010 

15.2 

22.5 

0 

NaN 

26.0 

38.1 

0 

NaN 

8.6 

10.8 

0 

NaN 

10172010 

144 

19.7 

0 

NaN 

25.1 

41.6 

0 

NaN 

8.9 

13.8 

0 

NaN 

11032010 

15.3 

27.0 

0 

NaN 

24.4 

36.9 

0 

NaN 

12.3 

16.0 

0 

NaN 

11172010 

20.2 

32.7 

0 

NaN 

30.2 

52.2 

0 

NaN 

14.1 

25.4 

3 

-0.9 

12032010 

15.9 

26.9 

0 

NaN 

25.6 

42.4 

0 

NaN 

11.6 

17.9 

a 

NaN 

12172010 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 


Table 21. Robustness of the Optimized Routes for San Diego to Guam 
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Tables 22-24 show what happens when the routes optimized with respect to a 
given environmental forecast encounter the best guess for the actual METOC conditions, 
which is the analysis environment. For Tables 22-24, routes are evaluated with respect 
to analysis, so the mean is given by: 


M' 


< ( 1 N 

^condition 

1=1 n =1 



j = analysis 


[ 10 ] 


For Tables 22-24, routes are evaluated with respect to analysis, so the max is given by: 
max | max \condition[w* n {i},j^u,j = analysis [11] 



Absolute Winds, kts(35ktThreshold) 

Relative Winds, kts (65kt Threshold) 

Seas, ft [25ft Threshold) 

Energy Consumption (khph) 


Departure Date 
IMMDDYYYY) 

mean 

max 

exceedances 

(number) 

exceedances 

Ikts) 

mean 

max 

exceedances 

(number) 

exceedances 

(kts) 

mean 

max 

exceedances 

(number) 

exceedances 

(ft) 

mean over 

route 

max over route 


0S012D10 

16.5 

25.5 

0 

NaN 

34,0 

43.1 

0 

NaN 

68 

9.8 

0 

NaN 

1841.8 

1849.9 


06152010 

18.4 

30.2 

0 

NaN 

32.2 

46.2 

0 

NaN 

5.5 

10.7 

0 

NaN 

1880.8 

1887.2 


07012010 

12.3 

19.6 

0 

NaN 

26.3 

38.8 

0 

NaN 

45 

7.2 

0 

NaN 

1822.9 

1827.4 

3 

07152010 

14.0 

23.5 

0 

NaN 

33.0 

47.8 

0 

NaN 

6.2 

8.6 

0 

NaN 

1842.5 

1846.2 

s. 

080I2DUO 

13.3 

23.7 

0 

NaN 

25.1 

42.1 

0 

NaN 

42 

8.1 

0 

NaN 

1819.2 

18218 

> 

08152010 

10.7 

16.3 

0 

NaN 

28.2 

44.5 

0 

NaN 

5.4 

7.4 

0 

NaN 

1853.5 

1860.7 

re 

c 

09012010 

13.0 

24.7 

0 

NaN 

316 

49.2 

0 

NaN 

as 

113 

0 

NaN 

1845.4 

1850.4 

■*-* 

09152010 

13.6 

22.0 

0 

NaN 

30.3 

50.1 

0 

NaN 

7.4 

149 

0 

NaN 

1884.0 

1900.5 

© 

F 

10012010 

14.5 

25.8 

0 

NaN 

29.2 

45.8 

0 

NaN 

9.2 

15.7 

0 

NaN 

1956.2 

1977.4 

c 

p 

10152010 

17.8 

29.6 

0 

NaN 

34.1 

54.2 

0 

NaN 

10.0 

147 

0 

NaN 

1902.3 

1919.8 

» 

11012010 

17.6 

30.0 

0 

NaN 

27.8 

42.1 

0 

NaN 

as 

142 

0 

NaN 

1940.9 

1984.8 


11152010 

12.7 

30.0 

0 

NaN 

23.7 

45.7 

0 

NaN 

as 

12.4 

0 

NaN 

1940.2 

1964.7 


12012010 

20.0 

30.7 

0 

NaN 

30.4 

53.6 

0 

NaN 

13.2 

25.9 

0.3 

** 

2328.2 

2466.7 


12152010 

18.1 

40.5 

0.4 

** 

30.6 

53.1 

0 

NaN 

9.8 

20.9 

0 

NaN 

1302.3 

2079.1 


Table 22. Optimized Routes for Norfolk to Rota Under Actual METOC Conditions 
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Absolute Winds, kts (35kt Threshold) 

Relative Winds, kts (65kt Thresh old) 

Seas, ft (25ftThreshold) 

Energy Consumption (khph) 


Departure Dale 
(MMDDYYYY) 

Mean 

Max 

I - 

*5 (o 

9 5 

< z 
m z: 

Exceedances 

(kts) 

c 

3 

2 

X 

ru 

2 

Exceedances 

(Number) 

8 

c 

re 

"O 

® 

8 2 
£ a 

Mean 

5 

2 

Exceedances 

(Number) 

8 

s 

■CT 

I 

8 — 
m ±i 
Lii 72 

Mean over 

Routs 

Maximum over 

Route 

j Environment-Analysis Route 

06022010 

14.1 

28.3 

0 

NaN 

28.1 

47.5 

0 

NaN 

5.8 

10.5 

0 

NaN 

2224.9 

2224.9 

(K162010 

** 

** 

** 

** 

** 

4r* 

** 

*4> 

*4< 

** 

** 

*4* 

*4< 

*4. 

07022010 

11.5 

24.5 

G 

NaN 

23.8 

41.7 

0 

NaN 

4.4 

8.6 

0 

NaN 

2079.4 

2224.1 

07162010 

*4r 

** 

** 

** 

** 

4** 

** 

4<4r 

4r4< 

** 

** 

** 

** 

*4* 

08022010 

** 

** 

** 

** 

** 

4** 

** 

*4» 

*4i 

** 

** 

*4* 

** 

*4* 

08162010 

*4> 

*4> 

#* 

** 

** 

4»* 

** 

*4r 

*4i 

** 

** 

*4r 

** 

*4> 

’09022010 

10.6 

22.2 

0 

NaN 

28.6 

48.0 

a 

NaN 

4.4 

9.6 

0 

NaN 

2198.2 

2223.7 

’09162010 

10.7 

20.7 

0 

NaN 

28.2 

40.8 

a 

NaN 

5.8 

8.1 

0 

NaN 

2268.5 

2275.7 

10022010 

*«r 

*4« 

*4t 

«* 

** 

4»* 

4<* 


*4< 

** 

** 

*4. 

** 

** 

10162010 

*4r 

** 

** 

** 

** 

** 

** 

*4» 

* * 


** 

*4r 

** 

** 

11022010 

15.5 

29.4 

c 

NaN 

247 

47.1 

0 

NaN 

9.9 

16.0 

0 

NaN 

2423.8 

2605.8 

’11162010 

15.0 

24.1 

0 

NaN 

29.8 

43.0 

0 

NaN 

10.8 

16.1 

0 

NaN 

2441.8 

2492.9 

'12022010 

** 

*4) 

** 

** 

** 

4r* 

** 

*4r 

4>4l 

** 

** 

** 

** 

*4> 

12162010 

** 

*4< 

** 

** 

*4< 

4>* 

** 

*4t 

4t 4t 

** 

** 

*4r 

** 

**> 


Table 23. Optimized Routes for Pearl Harbor to Yokosuka Under 

Actual METOC Conditions 




Absolute Winds, kts (35ktThreshold) 

Relative Winds, kts (65ktThresicld) 

Seas, ft (25ft Threshold) 

Energy Consumption (khph) 
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re 

2 

Exceedances 

(Number) 

8 

C 

re 

"O 

a 

s — 

x a 

LU 7L 

Mean oyer 

Route 

Maximum over 

Route 

Environment-Analysis Route 

06032010 

11.2 

20.3 

0 

NaN 

23.5 

41.1 

0 

NaN 

5.1 

10.8 

0 

NaN 

3149.6 

3599.6 

06172010 

** 

** 

**< 

** 

** 

»* 

»* 

** 

*4» 

** 

** 

** 

** 

4i4* 

07032010 

19.3 

32.7 

0 

NaN 

27.3 

47.8 

0 

NaN 

8.4 

13.8 

0 

NaN 

3729.5 

3729.5 

07172010 

Q 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Q 

0 

a 

08032010 

11.3 

16.4 

0 

NaN 

25.6 

38,3 

0 

NaN 

46 

6.8 

0 

NaN 

3538.8 

3551.5 

(38172010 

** 

** 

*4i 

** 

+* 

»* 

** 

M* 

*4» 

** 

4»4i 

*4t 

** 

4)4* 

09032010 

12.4 

24.1 

0 

NaN 

22.4 

49.6 

c 

NaN 

5.3 

9.4 

0 

NaN 

3373.2 

3636.3 

09172010 

14.6 

26.2 

0 

NaN 

24,8 

36.3 

0 

NaN 

7.9 

11.9 

0 

NaN 

3726.0 

3726.0 

10032010 

15.2 

31.2 

0 

NaN 

25.9 

55.1 

0 

NaN 

8.7 

12.6 

0 

NaN 

3781.7 

3846.6 

10172010 

16.7 

41.0 

1.3 

NaN 

26.2 

43.8 

0 

NaN 

8.5 

14.8 

0 

NaN 

3802.5 

3834.3 

11032010 

** 

4r* 

** 

** 

** 

>4 

** 

4»* 

4* 4* 

** 

4»* 

** 

** 

*4> 

11172010 

16.Q 

37.1 

1 

2.11 

26.3 

541 

0 

NaN 

11.1 

18.8 

0 

NaN 

3908.2 

3908.2 

12032010 

20.5 

40.5 

2 

5.40 

29.5 

57.4 

0 

NaN 

15.3 

31.5 

4 

1.4 

4180.0 

418C.Q 

12172010 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 


Table 24. Optimized Routes for San Diego to Guam Under 

Actual METOC Conditions 


Tables 25-27 show what happens when routes, optimized with respect to one forecast 
environment are evaluated with respect to environments represented by the 16 ensemble 
members. Throughout, means are taken with respect to all ensemble members, as in 
Equation 12: 


1 


M -1 



[ 12 ] 
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Similarly, the max values are maximum over the routes, as in Equation 13 with the 
optimization environment i given in the row, and max(stdev) taken over j - for 

all ensemble members. 

max ^ max [condition[w* n (z'),y'))j [13] 

In the top third of Tables 25-27, mean and standard deviation with respect to the 
environments for which the routes were optimized i are shown. 
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Absolute Winds, kts (35kt Threshold) 

Relative Winds, kts (65kt Threshold) 

Seas, ft (25ft Threshold) 

Energy Consumption 


Departure Date 
(MMDDYYYY) 

Mean 

Max 

Exceedances (Number) 

Exceedances(kts) 

Mean 

Max 

Exceedances (Number) 

Exceedances (kts) 

Mean 

Max 

Exceed a nces (N u m be r) 

Exceedances (ft) 

Mean over Route 

Maximum over Route 

All Ensemble Members Mean(StDev) ] 

06022010 

12.8(0.08) 

32.0(1.54) 

0.0(0.00) 

NaN 

27.8(0.23) 

49.2(0.47) 

0.0(0.00) 

NaN 

5.0(0.04) 

11.5(0.75) 

0.0(0.00) 

NaN 

2217.3(3.41) 

2256.4(3.62) 

06162010 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

07022010 

11.0(3.06) 

23.2(5.54) 

0.0(0.00) 

NaN 

24.3(6.54) 

42.7(10.52) 

0.0(0.00) 

NaN 

4.5(1.22) 

11.4(2.59) 

0.0(0.00) 

NaN 

2088.5(556.98) 

2286.2(563.34) 

07162010 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

08022010 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

08162010 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

09022010 

11.2(0.51) 

21.6(0.57) 

0.0(0.00) 

NaN 

28.4(0.42) 

54.8(4.07) 

0.0(0.00) 

NaN 

4.1(0.22) 

9.1(0.44) 

0.0(0.00) 

NaN 

2198.2(2.26) 

2235.0(5.73) 

09162010 

9.9(0.46) 

23.2(2.09) 

0.0(0.00) 

NaN 

27.5(0.34) 

43.1(0.05) 

0.0(0.00) 

NaN 

5.1(0.15) 

10.1(0.54) 

0.0(0.00) 

NaN 

2252.6(3.87) 

2285.2(5.43) 

10022010 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

10162010 

14.0(1.11) 

33.0(1.11) 

0.0(0.00) 

NaN 

29.3(0.19) 

51.6(0.90) 

0.0(0.00) 

NaN 

7.1(0.75) 

18.1(2.61) 

0.0(0.00) 

NaN 

2259.2(28.96) 

2340.0(45.11) 

11022010 

15.4(4.12) 

36.8(9.35) 

0.2(0.08) 

1.1(0.32) 

28.0(7.49) 

51.0(12.68) 

0.0(0.00) 

NaN 

10.2(2.79) 

17.4(4.31) 

0.0(0.00) 

NaN 

2443.2(654.06) 

2701.5(672.01) 

11162010 

16.4(1.28) 

43.5(3.67) 

0.8(0.22) 

4.4(0.28) 

30.2(0.82) 

66.7(6.09) 

0.1(0.02) 

1.7(0.00) 

11.9(0.70) 

27.7(2.91) 

0.1(0.05) 

2.7(0.00) 

2521.8(8.00) 

2662.0(11.77) 

12022010 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

12162010 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

*** 

| Analysis-Optimized Route j 

06022010 

12.6 

26.2 

0 

NaN 

28.4 

49.4 

0 

NaN 

4.8 

8.2 

0 

NaN 

2212.3 

2221.4 

06162010 

15.6 

31.7 

0 

NaN 

24.9 

40.2 

0 

NaN 

6.9 

14.1 

0 

NaN 

2324.1 

2385.7 

07022010 

11.5 

21.2 

0 

NaN 

26.6 

42.7 

0 

NaN 

4.7 

9.6 

0 

NaN 

2227.0 

2251.2 

07162010 

17.5 

90.0 

0.1 

55.0 

27.6 

67.6 

0.1 

2.6 

4.2 

8.4 

0 

NaN 

2198.2 

2219.7 

08022010 

12.5 

19.9 

0 

NaN 

29.9 

43.6 

0 

NaN 

4.6 

7.2 

0 

NaN 

2295.4 

2308.1 

08162010 

12.2 

24.5 

0 

NaN 

29.6 

43.1 

0 

NaN 

4.9 

11.0 

0 

NaN 

2262.7 

2289.6 

09022010 

11.1 

19.9 

0 

NaN 

28.6 

38.5 

0 

NaN 

4.0 

7.5 

0 

NaN 

2198.1 

2213.2 

09162010 

9.5 

18.8 

0 

NaN 

27.6 

43.1 

0 

NaN 

5.2 

10.1 

0 

NaN 

2254.2 

2269.5 

10022010 

13.6 

21.0 

0 

NaN 

27.6 

41.9 

0 

NaN 

5.8 

11.0 

0 

NaN 

2261.2 

2287.3 

10162010 

14.5 

30.6 

0 

NaN 

29.2 

51.6 

0 

NaN 

7.5 

18.1 

0 

NaN 

2273.2 

2340.0 

11022010 

16.5 

36.8 

0.2 

1.2 

29.6 

51.0 

0 

NaN 

11.1 

17.4 

0 

NaN 

2630.7 

2701.5 

11162010 

17.1 

43.5 

0.9 

4.7 

30.1 

66.7 

0.1 

1.7 

12.3 

27.7 

0.2 

2.7 

2503.3 

2615.0 

12022010 

24.7 

47.5 

1.9 

4.1 

39.1 

66.0 

0.1 

1.0 

13.9 

26.8 

0.7 

1.2 

2453.1 

2556.8 

12162010 

19.4 

44.8 

0.3 

3.3 

32.9 

63.1 

0 

NaN 

10.5 

24.9 

0 

NaN 

2398.9 

2466.8 

| Average-Environment-Optimized Route | 

06022010 

12.8 

32.4 

0 

NaN 

27.9 

49.4 

0 

NaN 

4.9 

11.7 

0 

NaN 

2215.4 

2257.7 

06162010 

15.6 

31.7 

0 

NaN 

24.9 

40.2 

0 

NaN 

6.9 

14.1 

0 

NaN 

2324.1 

2385.7 

07022010 

10.9 

22.0 

0 

NaN 

26.6 

42.7 

0 

NaN 

4.6 

9.7 

0 

NaN 

2221.9 

2238.0 

07162010 

12.9 

22.8 

0 

NaN 

28.3 

42.3 

0 

NaN 

4.3 

8.4 

0 

NaN 

2203.2 

2221.6 

08022010 

12.4 

18.8 

0 

NaN 

29.3 

45.9 

0 

NaN 

4.5 

7.1 

0 

NaN 

2224.7 

2235.5 

08162010 

12.3 

24.5 

0 

NaN 

29.6 

43.9 

0 

NaN 

4.9 

10.9 

0 

NaN 

2220.9 

2248.7 

09022010 

11.1 

19.9 

0 

NaN 

28.5 

38.5 

0 

NaN 

4.0 

7.5 

0 

NaN 

2198.1 

2212.8 

09162010 

9.8 

23.2 

0 

NaN 

27.4 

43.1 

0 

NaN 

5.0 

9.1 

0 

NaN 

2249.7 

2265.5 

10022010 

13.7 

25.0 

0 

NaN 

27.8 

41.9 

0 

NaN 

6.0 

12.1 

0 

NaN 

2270.6 

2297.2 

10162010 

14.5 

30.6 

0 

NaN 

29.2 

51.6 

0 

NaN 

7.5 

18.1 

0 

NaN 

2273.2 

2340.0 

11022010 

16.5 

36.8 

0.2 

1.2 

29.6 

51.0 

0 

NaN 

11.1 

17.4 

0 

NaN 

2630.7 

2701.5 

11162010 

16.5 

43.5 

0.6 

4.8 

30.5 

66.7 

0.1 

1.7 

11.9 

27.7 

0.2 

2.7 

2508.1 

2645.1 

12022010 

18.3 

35.7 

0.1 

0.7 

33.0 

55.8 

0 

NaN 

11.8 

26.0 

0.1 

1.0 

2509.9 

2666.0 

12162010 

19.3 

44.8 

0.3 

3.3 

32.9 

63.1 

0 

NaN 

10.5 

24.9 

0 

NaN 

2400.0 

2466.8 


Table 26. Overall Mean and Maximum Safety Experiences Along Pearl Harbor to 

Yokosuka Routes 


51 




































































nfc^ciuiDewrdSL kcs(35k: imeshcM) 

fldacve MMs. fesfcThrtshda] 

5«rc. "c(25ftnre»iicu) 

Ererx/Consu-nFOon ikfifli) 


a 

s 

s 

|! 

n 

E 

S 


1 

€ 

g 

I 

S 

V 

fi 

5 

j£ 

a 

s 

8 

S 


s 

6 

g 

a 

s 

XI 

8 

3 

M 

a 

s 

XI 

8 

s 


| 

I 

a 

5 

n 

f 

a 

i 

XI 

8 

9 

J 

i 

s 

a 

a 

o 

X 

m 

3 

| 

1 


a 

E 

i 

5 


a 

5 

S 

a 

a 

2 

2 

i 

a 

5 

5 


36C*2313 

1L 2(438 

22.7ULL2) 

3.3803) 


222(868] 

S41(154* 

13(3.33? 


6K23S] 

2L*7.27j 

00(13* 

HtH 

322L2( 1257.741 

3S611113169* 

S- 

36172313 

17.0(15* 

3LC( 7.721 

8.0(6.00) 


21*56?' 

490(12.24 

0.0(0331 

NeM 

8.3(115' 

14.5(153) 

00(10* 

HtH 

354L4t 94439) 

3ES2*970.1* 

Q 

ZTZ12Z.Z 

17.6(1081 

342(004) 

10(0.03) 

HtH 

263(004) 

47.8(10(1 

0.0(033) 

NsM 

0*005) 

158(0.16) 

00(10* 

HtH 

*7723(13() 

386L 4(871) 


vitaix 


2L4/.25) 

3.1ILUU) 

HtH 

Li 4866] 

39/03.54 

LUR1UU) 

Nffl 

i!*Lbl] 

/.Si*! 

U.SK1U* 

HtH 

3123.4(121110/) 

A/5.4I122L2* 

S 

seaizsio 

121(193 

26.4(287) 

1*600) 

HtH 

236(031) 

43.5(2.80 

00(033) 

HtH 

49(033) 

96(033) 

00(10* 

HtH 

35692(134) 

36546(107* 

s 

38172010 

8.2(218 

16.5(2. IS) 

10(3.03) 

HtH 

23 7(132! 

414(1300 

00(033) 

HtH 

15(093] 

8.9(2.11) 

00(10* 

HtH 

Si? 11'88137] 

3575.38362* 

k 


uorisa 

23.7(5.92) 

10(0.03) 

HtH 

222(594) 

496(12.531 

tom oo 

HtH 

SOUW) 

11*112' 

3.3*1321 

HtH 

34227*91273) 

36943923.24i 

! 

39i?2am 

13. S( 1081 

2v.SiO.19j 

lotcoa 

HtH 

23S(80S) 

42.5(113 

0.0(000 

HtH 

7.5(3.35] 

14.7(018) 

00(10* 

HtH 

*7105(149) 

37660(04* 

S 


16.2(173 

49.0(4.931 

HL23) 

MILOS] 

259iasu 

£-4(7.3* 

0.0(003) 

HtH 

8*036] 

212(124) 

00(10* 

HtH 

*787.6(22.94] 

3952.1(443* 

— 

mi? 2 om 

1S.SU i* 

34jgl S4) 

lotcaa 

HtH 

2S9(033) 

52.0(18* 

0.0(033) 

HtH 

9.1(3.15] 

17.0(0.12) 

00(10* 

HtH 

5836.3(2163] 

»9L4(4109 


110*2012 


*** 

»» 

»». 


*** 


... 

... 


»» 

.»• 

**• 


c 

11172010 

16.6(104 

37.5(057) 

IM 

L4(124| 

260(306) 

52.7(10* 

00(033) 

Me* 

lL*O04) 

208(03(9 

00(10* 

HcH 

3996.1(6.14) 

4225.5(35.0* 

< 

12017010 

7! TfUJI 

^=5(3 TO) 

i<nri7l 


nanim 

63 8(14* 

rntnn* 

HtH 

144(007) 

31 5(0 75) 


03(014) 

4176 0(5 97) 

4745 K? 73) 


12172313 


HtH 

Hit* 

HtH 

HtH 

N?fl 

HtH 

HtH 

ft?H 

HtH 

HtH 

HtH 

RtK 

N9M 


160*2012 

12J 

317 

3 

HtH 

254 

54.L 

3 

HtH 

7.3 

213 

: 

HtH 

3683.6 

38611 


361723m 

18.2 

313 

3 

HtH 

227 

49J 

3 

HtH 

8.6 

14.S 

: 

HtH 

3777.9 

3S827 


37312313 

176 

342 

3 

HtH 

36.3 

476 

3 

HtH 

8.6 

118 

: 

HtH 

37723 

396 14 

3 

37172310 

- 





- 

« 

•- 

** 

** 

- 

— 

— 

•• 


32312313 

12.7 

218 

3 

HcH 

222 

56.? 

0 

Mc« 

5.1 

11 

i 

HcH 

3566.1 

36323 

1 

[KT?70in 

S£ 

Iff. 

3 

HtH 

* 7 

W5 

3 


3 ? 

55 

2 

HtH 

3577 4 

3537 t 


090*2010 

143 

217 

3 

HtH 

23.6 

496 

3 


6.5 

118 

: 

HtH 

36402 

3694^ 

a. 

D91720m 

I3.S 

249 

3 

HtH 

23.5 

42.5 

3 


7.5 

14.7 

; 

HtH 

3710.7 

3«S Ji 

9 

13332313 

16.1 

46.6 

3.8 

(U 

36.2 

58.6 

3 


8.7 

3.3 

: 

HtH 

38114 

3987.5 

i 

13172313 

1S.9 

312 

ai 

11 

25.3 

48.1 

3 

rfeM 

10.1 

113 

: 

HtH 

3916.4 

4335.9 

i 

11312313 

16.5 

3L2 

3.2 

Lt 

34.1 

45.4 

3 

HtH 

12.9 

3.8 

; 

HtH 

40915 

4*209 


11172313 

166 

329 

33 

L4 

36.3 

527 

3 

HtH 

119 

218 

; 

HtH 

3906.1 

42255 


123*2313 

2L7 

56.9 

19 

4.1 

*3.6 

63.? 

3 

HtH 

14.4 

3LS 

16 

14 

41263 

4345.* 


121720m 

Ktn 


HtH 


ft?N 

HtH 

KtH 

HtH 

ztk 

HtH 

Nrfl 

HtH 

UtH 

HtH 


363*2313 

L2J 

317 

3 


25.5 

S4.L 

3 

HtH 

6.9 

2L3 

: 

HtH 

3677.3 

386LL 

2 

36172313 

18.2 

313 

3 


22.7 

492 

3 

HtH 

8.6 

14.5 

: 

HtH 

3777.9 

38827 

§ 

373*2313 

176 

342 

3 


36.3 

47J 

3 

HtH 

8.6 

118 

: 

HtH 

3772 3 

33614 

« 

37172313 

13.5 

212 

3 

HtH 

25.3 

39.7 

3 

HtH 

4.4 

7.9 

5 

HtH 

3566.1 

3575.4 

i 

37312313 

125 

218 

3 

HcH 

23.1 

43.5 

3 

HcH 

5.3 

11 

3 

HcH 

35626 

3624* 

£ 

381723m 

86 

16.6 

3 

HtH 

252 

39.5 

3 

HtH 

17 

IS 

: 

HtH 

3S27.4 

3537.* 

A 

393*2313 

14J 

217 

3 

HtH 

23.6 

49j6 

3 

HtH 

6.5 

118 

3 

HtH 

36402 

36942 

1 

39172313 

13.5 

£.9 

3 

HtH 

23. S 

42.5 

3 

HtH 

7.5 

14.7 

3 

HtH 

3713.7 

3^652 


133*2313 

16.* 

421 

34 

U 

36.4 

59.L 

3 

HtH 

9.3 

218 

: 

HtH 

37626 

3912? 

-s 

1UL.'2UW 

lb.4 

JLfa 

3 

HtH 

A. 11 

52J 

3 

HtH 

9.1 

112 

2 

HtH 

3HL99 

Stbb.* 


11312313 

16.5 

36.2 

32 

Li 

24.5 

45.4 

3 

HtH 

12.6 

213 

3 

HtH 

407L9 

41209 

? 

11172313 

166 

37.9 

33 

LA 

36.3 

52? 

3 

HtH 

LL9 

ZLS 

: 

HtH 

39961 

42255 

s 

123*2313 

2L7 

56.9 

19 

11 

*3.6 

63.? 

3 

HtH 

14.4 

315 

16 

14 

41263 

4345* 


12172013 

KfM 

NsH 

M«N 

HtH 


K?fl 

HtH 

HtH 

ffeM 

NeM 

HtH 

HtH 

UtH 

Ntfl 


Table 27. Overall Mean and Maximum Safety Experiences along San Diego to Guam 

Routes 


Although routes for most of the origin-destination pairs and departure dates are 
within the thresholds, there are some origin-destination pairs and departure dates with 
results that exceed the safety thresholds. Tables 19-21 shows exceedances experienced 
by all environmental conditions, absolute winds, relative winds and seas. For the most 
part, the exceedances were mild, five knots or less for winds and one foot or less for seas. 

Tables 22-24 show results within the safety thresholds. There were no 
exceedances of any of the environmental conditions. Tables 21, 22 and 27 have 
departure dates that are annotated by double asterisks (**). The results that were 
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generated for these departures dates were incomplete and removed from the table. The 
incompletion of results is to be expected because these departure dates lacked results in 
the tables in the sensitivity section. Table 14 shows the lack of data that is seen in Table 
23. Tables 16-18 show the lack of data that is seen in Table 24 and Table 27. 

For some routes and environments, STARS output produced errors in wind 
speeds. These errors can be seen in the results that were generated for seven departure 
dates in Tables 26 and one departure date in Table 27. Data for these routes and 
environments have been excluded for the respective departure dates and are annotated by 
triple asterisks (***). With the data errors removed, the remaining data show moderate 
exceedances in winds, less than 10 knots of winds, and extreme exceedances in sea state, 
up to 10 feet. 

From mild to extreme exceedances, the routes proved to be mostly robust when 
acted upon by environmental conditions. As to be expected, exceedances occurred more 
frequently in the late fall and early winter months. The extreme exceedances of nearly 10 
feet occurred in the month of December. The time of year and the weather that occurs 
during that time do have effects on ship routes. 

In summary, the rare and mild exceedances of safety thresholds indicate that 
routes are robust to uncertainty, as represented by an ensemble forecast. 
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IV. CONCLUSIONS 


A. SUMMARY 

This thesis set out to answer questions regarding the sensitivity and robustness of 
ship routes to environmental conditions and how important it is to take into consideration 
environmental uncertainty when determining a robust ship route. The sensitivity of the 
optimized routes showed that there were few routes that differed greatly in distance 
varying environmental forecasted conditions. The robustness of the optimized routes 
showed that the routes were safe under various environmental conditions. The time of 
year, with respect to seasons, is a major factor in the uncertainty of environmental 
conditions. The time of year and the forecasted weather conditions should definitely be 
taken into consideration when trying to determine an optimal ship route. 

Given that weather is a variable that cannot be controlled, it is important to 
understand how that uncertainty can affect a ship route. The optimized routes that were 
generated by STARS showed little sensitivity and were robust to the environmental 
conditions. Although the environmental conditions proved to have an impact on the 
optimized routes, the generated routes were, for the most part, safe. However, because 
forecasted weather is not a guarantee but an estimate of a possible occurrence, ship routes 
based on forecasted environmental conditions are subject to change with the changing 
conditions. The optimized routes can be used as a basis, but it should be clear that even 
an optimized route may need to be altered to maintain robustness. 

B. FUTURE WORK RECOMMENDATIONS 

With more time, a thesis that covers an entire calendar year will guarantee all 
seasons and all possible environmental conditions. Within the year of research, careful 
documentation of actual weather conditions will provide a reference for the generated 
routes. In addition to expanding the length of time, the use of an additional ship type 
could show if the routes acted the same on different types of ships of different tonnage 
and different fuel sources. 
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APPENDIX A 


SAMPLE OF STARS INPUT FILE 


REQUEST ID = SREQUESTID 
CLASSIFICATION = UNCLASSIFIED 
CAVEAT = NONE 

REQUEST TYPE = SREQUESTTYPE 
DESCRIPTION = San Diego CA to Guam 
PASSAGE = San Diego CA to Guam 


SHIP NAME = EFAS 
SHIP CLASS = DD963.DAT 
SHIP FOR DRAFT = 11.7 
SHIP AFT DRAFT = 12.0 
SHIP TRANSV GM = 2.6 
MAX HEAD SEA = 25.0 
MAX BEAM SEA = 25.0 
MAX FOLLOW SEA = 25.0 
MAX TRUE WIND = 35.0 
MAX REL WIND = 65.0 
MAX SPEED = 25.0 
MIN SPEED = 18.0 
MIN DIST 35 = 60 
MIN DIST 50 = 120 
DEPARTURE DATE = SDEPARTjUREDATE 
DEPARTURE TIME = SDEPARTURETIME 
ARRIVAL DATE = SARRIVALDATE 
ARRIVAL TIME = SARRIVALTIME 
WIND MODEL = NOGAPS 
WAVE MODEL = WW3_GLOBAL 
CURRENT MODEL = NCOM_GLOBAL 
UPPER BOUND WAYPOINTS = 4 
UB NUMBER = 01 
UB LATITUDE = 32.655575 
UB LONGITUDE = -117.237769 
UB NAV TYPE = GC 
UB NUMBER = 02 
UB LATITUDE = 51.224481 
UB LONGITUDE = -150.712078 
UB NAV TYPE = GC 
UB NUMBER = 03 
UB LATITUDE = 45.918669 
UB LONGITUDE = 158.866392 
UB NAV TYPE = GC 
UB NUMBER = 04 
UB LATITUDE = 13.442458 
UB LONGITUDE = 144.665872 
UB NAV TYPE = GC 
LOWER BOUND WAYPOINTS = 4 
LB NUMBER = 01 
LB LATITUDE = 32.655575 
LB LONGITUDE = -117.237769 
LB NAV TYPE = GC 
LB NUMBER = 02 
LB LATITUDE = 5.723406 
LB LONGITUDE = -135.341061 
LB NAV TYPE = GC 
LB NUMBER = 03 
LB LATITUDE = 1.097358 
LB LONGITUDE = 161.473958 
LB NAV TYPE = GC 
LB NUMBER = 04 
LB LATITUDE = 13.442458 
LB LONGITUDE = 144.665872 
LB NAV TYPE = GC 
END 



Environmental Data 
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APPENDIX B. SAMPLE OF STARS EXECUTIBLE FILE 


#!/b in/bash 

STARSDRIVER=/ITRI_AOTSR/otsr/stars_driver.2 

COLLECTRESULTS=/ITRI_AOTSR/otsr/zipSTARSKmlHtmlResults 

$STARSDRIVER \ 

Route SanDiego2Guam \ 

DepartureDate 20100603 \ 

ArrivalDate 20101231 \ 

RoutexEnvironment M A2 B2 C2 D2 E2 F2 Q1 R1 SI T1 U1 VI W1 XI Y1 Z1 analysis" 

\ 

WeaxEnvironment "A2 B2 C2 D2 E2 F2 Q1 R1 SI T1 U1 VI W1 XI Y1 Z1 analysis" \ 
OnCompletionNotify slhall@nps.edu 
$COLLECTRESULTS\ 

Route SanDiego2Guam \ 

DepartureDate 20100603 

exit 
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APPENDIX C. SAMPLE OF STARS HYPERTEXT MARKUP LANGUAGE OUTPUT FILE 


routex.20100601.NorfolkRota.A2 

REQUEST ID = EFAS 20100601 
REQUEST TYPE = ROUTEX 

DESCRIPTION = Naval Station Norfolk to US Naval Station Rota Spain 
PASSAGE = Naval Station Norfolk to US Naval Station Rota Spain 
SHIP NAME = EFAS 
DEPARTURE DATE = 06/01/2010 
DEPARTURE TIME = 12:00:00 


Point 

Dtg 

Lat 

Lon 

Ship 

Speed 

Ship 

Course 

Wind 

Speed 

Wind 

Direction 

SigWave 

Height 

Sea 

Height 

Sea 

Period 

Sea 

Direction 

Swell 

Height 

Swell 

Period 

Swell 

Direction 

Current 

Speed 

Current 

Direction 

Horse 

Power 

Distance 

1 

201006011200 

36.944 

-76.339 

25.00 

91.03 

17.07 

217.13 

2.76 

1.18 

4.28 

201.00 

2.49 

3.37 

200.00 

0.00 

0.00 

25.11 

35.26 

2 

201006011324 

36.931 

-75.604 

18.28 

131.01 

17.07 

217.13 

2.76 

1.18 

4.28 

201.00 

2.49 

3.37 

200.00 

0.43 

63.43 

19.63 

54.26 

3 

201006011622 

36.334 

-74.757 

19.24 

56.39 

19.22 

208.38 

4.68 

3.87 

4.97 

202.00 

2.62 

4.07 

196.00 

0.61 

18.43 

11.54 

31.20 

4 

201006011800 

36.621 

-74.217 

19.24 

56.71 

20.49 

213.39 

6.01 

5.12 

5.34 

210.00 

3.15 

4.56 

212.00 

0.19 

90.00 

2.91 

7.16 

5 

201006011822 

36.686 

-74.093 

19.45 

136.71 

20.49 

213.39 

6.01 

5.12 

5.34 

210.00 

3.15 

4.56 

212.00 

0.19 

90.00 

27.28 

61.47 

6 

201006012131 

35.938 

-73.225 

18.78 

59.53 

19.17 

210.47 

5.57 

4.33 

5.12 

214.00 

3.51 

4.96 

220.00 

2.75 

45.00 

9.17 

37.73 

7 

201006012332 

36.254 

-72.553 

18.83 

53.32 

19.17 

210.47 

5.57 

4.33 

5.12 

214.00 

3.51 

4.96 

220.00 

2.92 

53.13 

2.06 

8.64 

8 

201006020000 

36.340 

-72.410 

18.83 

53.40 

19.64 

211.64 

5.60 

4.46 

5.18 

219.00 

3.38 

4.88 

230.00 

2.92 

53.13 

7.41 

31.08 

9 

201006020139 

36.648 

-71.892 

20.21 

43.99 

20.68 

220.43 

6.66 

5.41 

5.46 

223.00 

3.87 

5.39 

231.00 

2.44 

61.39 

13.91 

44.22 

10 

201006020350 

37.177 

-71.249 

18.74 

44.93 

20.50 

227.69 

7.35 

5.58 

5.57 

227.00 

4.79 

6.02 

240.00 

2.15 

84.81 

12.05 

40.50 
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11 

201006020600 

37.653 

- 70.647 

18.74 

45.31 

21.27 

233.55 

8.40 

[?.77 

5.53 

230.00 

6.10 

6.70 

237.00 

[ 2.17 

100.30 

11.03 

3.17 

12 

201006020610 

( 37.690 

- 70.600 

19.61 

34.25 

21.27 

' 233.55 

8.40 

5.77 

5.53 

230.00 

6.10 

6.70 

237.00 

2.17 

100.30 

17.95 

51.93 

13 

201006020849 

1 : 38.404 

- 69.978 

25.00 

81.90 

20.34 

231.60 

8.40 

5.28 

5.47 

230.00 

6.53 

6.70 

242.00 

0.19 

270.00 

25.88 

34.55 

14 

201006021011 

38.483 

- 69.250 

25.00 

81.79 

19.09 

234.53 

8.09 

4.07 

5.02 

232.00 

6.99 

6.69 

245.00 

0.78 

0.00 

23.22 

34.57 

15 

201006021134 

38.563 

- 68.520 

25.00 

65.99 

14.84 

248.48 

8.15 

2.56 

4.55 

246.00 

7.74 

6.62 

239.00 

0.58 

0.00 

7.18 

10.45 

16 

201006021200 

38.633 

- 68.317 

25.00 

66.12 

15.95 

247.81 

8.37 

2.92 

4.61 

245.00 

7.84 

6.74 

244.00 

0.58 

0.00 

17.69 

25.69 

17 

201006021301 

38.806 

- 67.814 

25.00 

81.94 

15.95 

247.81 

8.37 

2.92 

4.61 

245.00 

7.84 

6.74 

244.00 

0.19 

0.00 

24.91 

34.55 

| l 8 

201006021424 

38.884 

- 67.082 

25.00 

81.85 

17.62 

247.96 

8.64 

( 3.77 

4.94 

245.00 

7.78 

6.85 

249.00 

0.27 

225.00 

26.33 

34.56 

19 

201006021547 

38.964 

- 66.349 

18.88 

118.40 

16.78 

243.14 

8.32 

3.22 

4.75 

241.00 

7.68 

6.78 

250.00 

0.87 

243.43 

16.87 

41.68 

20 

201006021800 

38.631 

- 65.566 

18.88 

118.88 

16.78 

243.14 

8.32 

3.22 

4.75 

241.00 

7.68 

6.78 

250.00 

0.78 

0.00 

0.16 

0.43 

21 

201006021801 

38.627 

- 65.558 

25.00 

81.57 

16.78 

243.14 

8.32 

3.22 

4.75 

241.00 

7.68 

6.78 

250.00 

0.78 

0.00 

23.23 

34.70 

22 

201006021924 

38.710 

- 64.825 

25.00 

67.18 

17.82 

243.43 

8.64 

[ 3.84 

4.92 

243.00 

7.74 

6.95 

254.00 

1.13 

30.96 

23.28 

36.24 

23 

201006022051 

38.942 

- 64.109 

25.00 

81.64 

20.11 

249.64 

9.16 

4.89 

5.32 

248.00 

7.74 

7.14 

259.00 

0.80 

104.04 

26.08 

34.70 

24 

201006022214 

39.024 

- 63.373 

25.00 

81.62 

21.50 

255.34 

9.30 

5.05 

5.37 

253.00 

7.81 

[ 7.36 

261.00 

1.17 

0.00 

22.27 

34.72 

25 

201006022338 

39.106 

- 62.635 

18.39 

118.46 

21.50 

255.34 

9.30 

5.05 

5.37 

253.00 

7.81 

| 7.36 

261.00 

0.87 

( 333.43 

2.23 

6.67 

26 

201006030000 

39.052 

- 62.509 

18.39 

118.54 

21.50 

255.34 

9.30 

5.05 

5.37 

253.00 

7.81 

7.36 

261.00 

0.87 

333.43 

11.71 

35.08 

27 

201006030154 

38.771 

- 61.850 

18.69 

123.13 

20.13 

254.89 

9.93 

5.02 

5.35 

254.00 

8.56 

7.50 

266.00 

0.58 

0.00 

16.60 

44.34 

28 

201006030416 

38.365 

- 61.061 

25.00 

76.44 

20.11 

246.66 

9.17 

4.76 

5.29 

247.00 

7.84 

7.60 

277.00 

0.61 

198.43 

17.18 

21.26 

[29 

201006030507 

38.447 

- 60.621 

25.00 

100.34 

20.11 

246.66 

9.17 

4.76 

5.29 

247.00 

7.84 

7.60 

277.00 

0.39 

0.00 

15.40 

21.74 

30 

201006030600 

( 38.381 

- 60.167 

25.00 

100.62 

19.93 

249.44 

9.32 

4.66 

5.27 

249.00 

8.07 

( 7.77 

283.00 

0.80 

14.04 

10.21 

15.06 

I 31 

201006030636 

38.335 

- 59.852 

18.20 

52.64 

19.93 

249.44 

9.32 

4.66 

5.27 

249.00 

8.07 

7.77 

283.00 

0.80 

14.04 

15.14 

44.15 

32 

201006030901 

( 38.779 

- 59.102 

25.00 

80.54 

15.61 

251.11 

8.95 

2.82 

4.62 

249.00 

8.50 

7.53 

276.00 

0.99 

101.31 

27.09 

36.33 

33 

201006031028 

38.876 

- 58.335 

18.44 

48.24 

16.08 

255.29 

9.05 

2.92 

4.65 

253.00 

8.56 

7.67 

281.00 

[ T .05 

111.80 

10.49 

28.00 

34 

201006031200 

39.186 

- 57.886 

18.44 

48.53 

16.08 

255.29 

9.05 

2.92 

4.65 

253.00 

8.56 

7.67 

281.00 

0.97 

126.87 

7.34 

18.70 

[35 

201006031300 

39.391 

- 57.583 

25.00 

87.33 

16.08 

255.29 

9.05 

2.92 

4.65 

253.00 

8.56 

7.67 

281.00 

0.87 

[ 63.43 

25.32 

36.02 

36 

201006031427 

39.417 

- 56.807 

18.15 

48.16 

16.64 

255.80 

9.18 

2.95 

4.66 

253.00 

8.69 

7.73 

285.00 

0.97 

36.87 

15.29 

46.66 

I 37 

201006031701 

39.933 

- 56.051 

18.21 

48.20 

13.16 

257.20 

8.63 

[1.31 

4.14 

253.00 

8.53 

7.70 

291.00 

11.23 

71.57 

5.69 

17.74 

[38 

201006031800 

40.130 

- 55.763 

18.21 

48.39 

13.16 

257.20 

8.63 

1.31 

4.14 

253.00 

8.53 

7.70 

291.00 

1.23 

71.57 

9.20 

28.68 

39 

201006031934 

40.446 

- 55.294 

18.65 

122.16 

12.73 

262.98 

8.93 

1.15 

4.14 

261.00 

8.86 

7.67 

294.00 

0.61 

18.43 

15.88 

43.76 

40 

201006032155 

40.055 

- 54.487 

18.39 

48.19 

13.77 

256.94 

8.39 

1.61 

4.31 

253.00 

8.23 

7.63 

298.00 

0.97 

716.87 

12.80 

38.23 

41 

201006040000 

40.479 

- 53.862 

18.39 

48.60 

13.77 

256.94 

8.39 

1.61 

4.31 

253.00 

8.23 

7.63 

298.00 

0.97 

53.13 

2.83 

8.44 


62 
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